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SECTION  I  1 

INTRODUCTION  I 


The  purpose  of  the  research  presented  in  this  document  is  to  refine  our  ability  to 
discriminate  between  nuclear  explosions  and  earthquakes  based  on  the  depth  of  an  event 
determined  from  local  and  regional  phases.  If  the  depths  of  seismic  events  can  be  deter¬ 
mined  accurately,  many  events  can  be  eliminated  as  potential  nuclear  explosions  based 
on  the  practical  limit  of  nuclear  explosion  burial  depth.  This  discriminant  has  been  used 
for  many  years  at  teleseismic  distances.  The  higher  frequency  data  obtained  at  local  and 
regional  distances  improves  the  depth  resolution  of  the  method  and  thus  the  minimum 
depth  at  which  the  discriminant  may  be  applied.  In  addition,  smaller  events  not 
recorded  at  teleseismic  distances  may  be  seen  at  regional  and/or  local  distances. 

Our  approach  to  this  problem  has  been  to  improve  phase  identification  and  picking 
abilities  of  depth  phases  pP  and  aP  through  filtering  and  three-component  polarization 
methods.  These  methods  are  applied  first  to  data  from  a  set  of  earthquakes  located  in 
the  northeastern  United  States  and  adjacent  Canada.  The  area  of  study  was  chosen  for 
several  reasons: 

(1)  the  large  majority  of  events  in  this  area  occur  at  depths  less  than  15  km  so  that  a 
data  set  of  shallow  events  is  readily  available, 

(2)  the  geology  and  tectonic  setting  of  the  area  is  well  known  and  similar  to  that  found 
in  the  Soviet  Union, 

(3)  the  depths  of  the  events  in  the  area  of  study  are  fairly  well  known  from  network 
studies  and  aftershock  surveys,  and 

(4)  there  are  two  broad-band  three-component  digital  seismic  stations  in  the  vicinity: 
RSNY  and  SRNY.  RSNY,  part  of  the  Regional  Seismic  Test  Network  (RSTN),  is 
in  northern  New  York  and  SRNY,  operated  by  Rondout  Associates,  is  in 
southeastern  New  York. 

Having  refined  the  techniques  for  depth  determination  in  the  well  known  region  of 
the  northeastern  U.S.,  the  methods  were  applied  to  a  set  of  data  obtained  from  the 
Ocean  Sub-bottom  Seismometer  (OSS)  installed  off  the  Kuril  trench  in  the  north  Pacific 
Ocean.  Several  events  from  the  Kamchatka  peninsula  and  the  Kuril  Islands  were  stu¬ 
died. 

The  most  difficult  problem  associated  with  depth  determination  at  local  and 
regional  distances  using  depth  phases  is  identification.  As  the  depth  phases  are  usually 
buried  in  the  coda,  data  processing  plays  a  key  role  in  identifying  and  picking  them. 
Band-pass  filtering  and  three-component  polarization  state  filtering  methods  are  used  on 
the  data  in  an  attempt  to  isolate  the  depth  phases  and  adaptive  polarization  analysis  is 
used  to  help  identify  them. 

Once  the  depth  phases  have  been  correctly  identified  and  picked,  the  depth  can  be 
determined  accurately  if  the  source  region  velocity  structure  is  known.  However, 
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misidentification  of  phases  can  lead  to  severe  errors  in  the  depth  estimates.  To  assess  the 
potential  depth  errors  resulting  from  these  problems,  synthetic  seismogram  depth  sec¬ 
tions  have  been  generated  for  each  of  the  events  in  the  northeastern  U.S.  to  compare  to 
the  identification  and  time  picks  made  from  the  data.  Of  course,  the  depth  of  the  event 
is  not  the  only  factor  contributing  to  the  character  of  the  depth  section;  the  velocity- 
depth  structure  used  and  the  focal  mechanism  are  also  very  important  to  the  overall 
character  of  the  section. 

In  the  sections  that  follow,  we  briefly  describe  the  analysis  methods  used  on  the 
data  to  isolate  depth  phases,  review  the  set  of  earthquakes  that  we  have  selected  for 
study,  and  show  detailed  analyses  of  the  earthquakes. 
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l 


EVENT  DEPTH  DETERMINATION 


In  this  report  we  describe  a  method  of  depth  determination  that  involves  enhancing 
the  suspected  depth  phases  at  the  beginning  of  the  P  coda  using  advanced  polarization 
methods  on  three-component  broad-band  data.  Single  station  techniques  were 
emphasised  as,  in  a  realistic  test  ban  monitoring  situation,  it  is  likely  that  only  one  sta¬ 
tion  will  record  many  of  the  suspect  events.  The  particular  techniques  used  to  enhance 
the  depth  phases  were  adaptive  polarization  analysis,  a  method  that  adaptively  rotates 
the  horizontal  components  to  the  direction  of  maximum  and  minimum  energy  propaga¬ 
tion,  and  polarization  state  filtering  (Samson  and  Olsen,  1081),  an  adaptive  frequency 
filtering  technique  that  passes  only  those  frequencies  in  a  specified  polarisation  state  (e.g. 
rectilinear).  These  techniques  successfully  enhance  polarised  energy  relative  to  the  scat¬ 
tered  component  of  the  coda,  which  typically  complicates  the  interpretation  of  local  and 
regional  seismograms.  We  note,  however,  that  our  research  has  concentrated  primarily 
on  events  with  signal- to-noise  ratios  of  one  or  better,  and  studies  on  the  effectiveness  of 
the  polarization  state  filter  have  shown  that  its  effectiveness  deteriorates  below  this  level. 
Therefore,  if  depth  phases  are  to  be  used  to  determine  depth  on  very  small  magnitude 
events,  then  either  better  methods  of  enhancing  the  depth  phases  must  be  found  or 
receiver  sites  must  have  low  noise  characteristics. 

The  next  phase  of  the  processing  scheme  is  to  identify  the  individual  depth  phases 
either  through  experience,  synthetic  seismogram  comparison,  or  some  other  method. 
Because  the  depth  phases  all  have  the  same  phase  velocity  and  polarisation  characteris¬ 
tics,  identification  of  the  phases  cannot  be  made  on  these  parameters  alone.  Additional 
information  such  as  relative  arrival  times  and  amplitudes  must  be  used.  Whatever 
scheme  is  used  to  identify  the  arrivals,  it  is  generally  assumed  that  the  source  is  a  simple 
point  mechanism  within  the  frequency  band  being  studied.  For  events  at  local  and 
regional  distances,  this  assumption  is  not  always  valid. 

Polarisation  State  Filtering 

Polarisation  analysis  of  a  three-component  wavefield  is  used  to  determine  the  vector 
particle  motion  characteristics  such  as  the  degree  of  polarisation,  direction  of  propaga¬ 
tion  and  ellipticity.  This  information  can  then  be  used  to  design  filters  that  emphasize 
or  eliminate  waves  having  specific  polarization  properties.  In  seismology,  this  method 
has  been  applied  to  analyse  and  filter  waves  coming  from  certain  directions,  or  having 
rectilinear  or  elliptical  particle  motion,  and  to  discriminate  against  isotropic  noise  (Sut¬ 
ton  and  Pomeroy,  1983;  White,  1984;  Archambeau  et  al. ,  1985)  . 

Polarized  arrivals  can  be  enhanced,  their  wave  type  identified,  and  their  direction  of 
propagation,  as  well  as  their  phase  velocity  for  body  waves  and  ellipticity  for  Rayleigh 
waves,  determined.  This  information  can  then  be  used  for  phase  identification,  epicen¬ 
tre!  location,  source  depth  determination,  earth  structure  studies,  ...etc.  Data  filtered  for 
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polarised  signal  are  more  amenable  to  comparison  with  synthetic  seismograms,  providing 
a  better  determination  of  relevant  parameters. 

The  polarisation  state  method  was  proposed  by  Samson  (1077)  for  the  analysis  and 
filtering  of  seismograms.  It  makes  use  of  the  spectral  matrix,  which  is  the  frequency 
domain  equivalent  of  the  cross-correlation  matrix  of  the  data,  smoothed  over  a  frequency 
window.  Various  spectral  estimators  can  be  obtained  from  the  spectral  matrix,  such  as 
the  degree  of  polarisation,  the  degree  of  rectilinear  polarisation,  a  detector  for  signal  in  a 
specific  "pure  state"  (i.e.  completely  polarised)  or  with  specific  polarisation  characteris¬ 
tics.  Such  estimators  can  then  be  used  to  filter  the  data  in  order  to  discriminate  against 
unpolarised  noise  or  to  enhance  signal  in  some  specific  polarisation  state  (Samson  and 
Olson,  1081).  Since  the  polarisation  characteristics  of  seismograms  continuously  change 
with  time,  such  filters  have  to  be  data- adaptive:  a  sliding  time  window  of  data  is 
analysed  and  a  filter  is  specifically  designed  for  that  window.  Polarisation  properties  are 
determined  and  used  not  only  as  a  function  of  time,  but  also  as  a  function  of  frequency, 
since  the  processing  is  done  in  the  frequency  domain.  This  is  an  advantage  that  methods 
which  operate  purely  in  the  time  domain  do  not  have. 

To  enhance  depth  phases  pP  and  »P  we  take  advantage  of  the  three-component 
polarisation  state  filter’s  ability  to  pass  frequencies  that  contain  rectilinearly  polarised 
motion  and  reject  frequencies  that  do  not. 

After  the  data  have  been  state  filtered  for  rectilinear  motion,  they  are  processed 
with  an  adaptive  polarisation  method.  The  horisontal  components  are  adaptively 
rotated  to  the  time  varying  radial  and  transverse  directions  and  the  time  variable 
azimuth  and  apparent  angle  of  incidence  are  given.  This  process  helps  in  the 
identification  of  phases.  Also,  the  synthetic  seismogram  components  are  the  radial  and 
vertical  so  rotating  the  horisontal  data  to  the  adaptive  radial  and  transverse  components 
facilitates  comparison. 

Adaptive  Polarisation  Analysis 

In  the  Adaptive  Polarisation  Analysis  method  developed  at  RAI,  the  horisontal 
components  are  adaptively  rotated  to  the  time  varying  radial  and  transverse  directions 
and  the  time  variable  asimuth  and  apparent  angle  of  incidence  are  given.  This  process 
helps  in  the  identification  of  phases.  Below  we  present  the  theory  behind  this  method 
assuming  rectilinear  motion  at  a  constant  asimuth  and  angle  of  incidence.  Later  in  the 
discussion,  we  will  consider  simultaneous  arrivals  from  different  directions. 

We  seek  to  determine  the  asimuth  and  angle  of  incidence  over  a  small  time  window; 
one  that  is  at  least  as  long  as  one  or  two  cycles  of  the  predominant  frequency.  The 
asimuth  is  determined  by  examining  the  sero  lag  cross  covariance  between  the  two  hor¬ 
isontal  components  as  a  function  of  azimuth  angle.  One  of  the  angles  at  which  this  func¬ 
tion  is  sero  is  the  asimuth  of  the  predominant  particle  motion.  The  horisontal  com¬ 
ponents  of  motion  as  functions  of  the  north  (N)  and  east  (E)  horizontals  and  the 
asimuth  (9)  are  written: 


■ 
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Hj(t)  -  N(t)eos0  +  E(t)sin0 

Hj(t)  -  E(t)cos0  -  N(t)ain0  (1) 

and,  the  cross-covariance  between  Hj  and  H,  at  sero  lag  is: 

1  ■ 

ccv0  *  S  S,(t+iAt)  Sj(t+iAt)  (2) 

n  i-o 

where  fi(t)  »  H(t)  ~  Mg,  the  mean  value  of  H(t)  over  the  interval,  and  At  is  the 
sampling  interval.  In  terms  of  17  and  £: 


ecvo(0)  “  a  cos  20  +  6  sin  20  (3) 

where 

1  » 

a  —  — SW(t  +  iAt)E(t  +  iAt) 

“  i-o 

1  *  ,  , 
ft  -  —  S(E  (t  +  iAt)  -  N  (t  +  iAt)] / 2 

“i-o 

Solving  for  0  at  ccvo0  »  0: 

0  -  tan-l(-a/b)/2  (4) 

• 

The  azimuth  0  obtained  from  (4)  will  be  between  -45  *  and  45  *  and  could  represent 
any  one  of  four  directions;  the  direction  of  particle  motion,  its  supplement,  its  transverse, 
or  the  transverse  supplement.  Without  using  the  vertical  component,  the  direction 
obtained  from  (4)  can  be  determined  with  a  180  *  ambiguity  by  examining  the  behavior 
of  (3).  The  sero  lag  croaa-covariance  as  a  function  of  azimuth  given  by  (3)  is  a  sinusoid. 
Suppose  that  the  0  determined  from  (4)  aligns  the  H,  component  with  the  direction  of 
maximum  particle  motion.  Increasing  this  #  in  (3)  will  place  the  direction  of  particle 
motion  between  the  Hj  and  H,  components  and  the  value  of  ccv0  will  increase.  Con¬ 
versely,  if  Hj  is  aligned  with  the  direction  of  particle  motion  by  (4)  then  increasing  0  in 
(3)  will  decrease  ccvQ.  Thus,  when  the  slope  of  the  sinusoid  given  by  (3)  at  the  azimuth 
determined  in  (4)  is  positive,  then  f  is  transverse  to  the  direction  of  particle  motion  (or 
its  supplement)  and  when  the  slope  of  the  sinusoid  given  by  (3)  at  the  azimuth  deter¬ 
mined  in  (4)  is  negative,  then  0  is  the  direction  of  particle  motion  (or  its  supplement). 
Rather  than  computing  the  slope  from  the  derivative  of  (3)  with  respect  to  0,  we  can  find 
its  sign  by  examining  (3)  at  0  «  0  where  ccvQ(0)  =  a,  If  the  signs  of  a  and  the  0  found 
from  (4)  are  similar  (i.e.  both  have  the  same  sign),  then  the  slope  is  negative  and  0  is  the 
direction  of  particle  motion  or  its  supplement.  If  the  signs  are  dissimilar,  then  the  slope 
is  positive  and  the  direction  of  particle  motion  is  at  either  0  +  90 '  or  0  -  90  * .  The 
direction  of  particle  motion  can  be  either  longitudinal  or  transverse  to  the  propagation 
direction  depending  on  whether  the  particle  motion  is  P/SV  or  SH.  If  the  azimuth  from 
the  event  to  the  station  is  known,  then  particle  motion  between  ±  45  *  of  that  azimuth 
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can  be  considered  radial  and  between  45  *  and  135  *  or  —45  *  and  —135  *  can  be  con¬ 
sidered  transverse.  The  asimuth  obtained  from  (4)  together  with  the  assumed  azimuth  to 
the  event  can  be  used  to  rotate  the  horizontal  components  to  the  radial  and  transverse 
directions  at  the  center  of  the  window.  The  window  is  then  advanced  one  sample  and 
the  process  repeated.  The  radial  and  transverse  traces  separate  the  P/SV/Rayleigh 
motion  from  the  SHf  Love  motion. 

The  apparent  angle  of  incidence  is  calculated  in  much  the  same  manner  as  the 
asimuth.  For  this  measurement,  the  adapted  radial  and  vertical  components  are  used  as 
the  orthogonal  components  from  which  the  angle  of  incidence  is  determined.  The  dis¬ 
tinctive  particle  motions  of  P,  SV,  and  Rayleigh  waves  are  distinguished  using  the 
apparent  angle  of  incidence.  P  motion  gives  a  positive  value,  SV  a  negative  one  and 
Rayleigh  a  rotating  angle  of  incidence. 

An  estimate  of  the  relative  accuracy  of  the  asimuth  and  angle  of  incidence  measure¬ 
ments  is  obtained  through  the  maximum  value  of  the  zero  lag  cross-correlation  function. 
At  its  maximum  value,  where  the  predominant  direction  of  motion  bisects  Ht  and  H2, 
the  zero  lag  cross-correlation  function  (cef0(#  +  45  * ))  will  be  between  0  and  1,  inclusive. 
A  large  value  is  indicative  of  strong  rectilinear  motion  over  the  window  and  a  small 
value  indicates  the  lack  thereof.  Where  there  is  little  or  no  rectilinear  motion,  confidence 
in  the  results  obtained  is  also  small  and  we  expect  the  error  to  be  close  to  00  * .  Con¬ 
versely,  a  large  value  for  ccf0(f  +  45  * )  indicates  strong  rectilinear  motion  and  should 
have  an  error  close  to  0  * .  This  subjective  error  evaluation  can  be  written: 

error  *  90  *  (1  -  ccf0(f  +  45  * )).  (5) 

Error  determined  in  this  manner  is  not  an  absolute  quantity  but  rather  a  relative 
measure  of  the  accuracy  of  the  value. 

Because  the  method  presented  here  relies  on  the  determination  of  the  predominant 
direction  of  motion,  the  question  arises;  "What  angle  would  be  determined  if  two 
uncorrelated  signals  arrived  at  the  station  at  the  same  time  from  different  directions?". 
Let  the  two  uncorrelated  signals  be  x(t)  and  y(t)  where  the  standard  deviations  of  x  and 
y  are  equal  and  define  the  north  and  east  components  of  motion  as 

N(t)  -  jx(t)  +  ky(t) 

E(t)  -  rx(t)  +  sy(t)  1 

j,  k,  r,  and  s  are  scaling  factors  that  can  be  set  to  simulate  any  combination  of 
direction  and  amplitude  for  x  and  y.  Substituting  the  values  for  N  and  E  in  (6)  into 
equation  (4) 
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ElO*  +  ky)(rx  +  sy)] 

S((rV  +  «V  +  2my- jV  -  kV  -  2jkxy)]/2 


(T) 


Because  x  and  y  are  not  correlated,  the  xy  terms  tend  to  sero  as  the  length  of  the 
series  gets  large.  Rewriting  (7)  eliminating  these  terms 

1  . 

9  *  — tan  • 

2 

Equation  (8)  may  be  used  to  predict  9  for  two  time  series  passing  a  station  at 
different  asimuths.  The  worst  case  occurs  when  the  two  wavetrains  are  traveling  at 
asimuths  45  *  apart,  for  when  one  of  the  signals  is  showing  no  correlation  at  its  asimuth, 
the  other  will  be  at  a  maximum.  If  the  two  signals  are  orthogonal,  then  the  proper 
asimuth  would  be  chosen. 

The  above  analysis  is  valid  only  in  the  case  that  the  two  signals  are  not  coherent. 
When  the  two  signals  are  perfectly  coherent  then  the  asimuth  determined  will  be  that  of 
the  vector  sum  of  the  two  signals.  Partial  coherency  would  result  in  some  combination 
of  the  completely  non-coherent  result  and  the  perfectly  coherent  result. 

In  Figure  1,  we  show  synthetic  vertical  and  horisontal  components  to  be  used  as  the 
input  to  the  polarisation  program.  The  modulated  arrivals  on  the  horisontal  com¬ 
ponents  represent  a  constant  amplitude  signal  that  rotates  smoothly  in  asimuth.  A  con¬ 
stant  amplitude  vertical  component  maintains  an  unchanging  angle  of  incidence.  The 
results  of  processing  this  data  with  the  polarisation  program,  “trvming  that  the  signal 
contains  longitudinal  motion  exclusively,  are  shown  in  Figure  2A.  There  are  six  output 
traces;  the  vertical,  adapted  radial,  and  adapted  transverse  components  of  motion,  the 
product  of  the  radial  and  vertical  traces,  and  the  asimuth  and  apparent  angle  of 
incidence  traces  plotted  with  error  bars.  The  radial  trace  exhibits  a  constant  amplitude 
and  the  transverse  component  shows  very  little  motion  at  all.  Thus,  the  radial  follows 
the  predominant  direction  of  motion.  In  the  case  of  a  short  duration  arrival,  the  RZ 
trace  would  help  to  identify  the  first  arrival  and  its  direction  but  in  this  ease,  where  the 
signal  is  continuous,  RZ  has  little  use.  The  asimuth  trace  tracks  the  motion  of  the  rotat¬ 
ing  arrival  quite  well  and  shows  very  little  error.  The  same  is  true  of  the  apparent  angle 
of  incidence  trace  where  there  is  no  ehange  throughout  the  duration  of  the  analysis. 
When  a  propagation  direction  is  sssumed  in  the  polarisation  program  (Figure  2b),  the 
particle  motion  jumps  between  the  radial  and  transverse  components  ss  the  asimuth 
rotates  from  one  direction  to  the  other.  The  RZ  and  apparent  angle  of  incidence  traces 
help  to  identify  the  type  of  motion.  When  positive,  the  motion  is  P  type;  when  nega¬ 
tive,  SV\  and  when  sero,  SH.  Rayleigh  waves  produce  a  rotation  in  the  RZ  and  angle  of 
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Vertical,  North-South,  and  East-West  com¬ 
ponents  of  motion  for  a  signal  sweeping  3G0  *  of 
asimuth  plus  noise. 
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Figure  tA.  Adaptive  polarisation  analysis  of  the  data  in 
Figure  1  Morning  no  asimuth.  Z— vertical; 
R— radial;  T— transverse;  RZ— sign  (RZ)  \/RZ; 
AZ— asimuth  with  relative  error  bars;  AA— angle 
of  incidence  with  error  bars.  The  radial  trace 
tracks  the  direction  of  particle  motion  around 
the  circle. 


seconds 


Figure  tB.  Same  u  Figure  2A.  but  assuming  an  asimuth 
0  * .  The  motion  changes  from  radial 
transverse  and  back  as  the  azimuth  moves 
around  the  circle.  Positive  values  of  RZ  and  AA 
indicate  P  type  motion;  zero  values  indicate  SH 
motion;  and  negative  values,  SV  motion. 
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incidence  traces.  Note  that  the  asimuth  trace  jumps  from  +45*  to  -45*  as  the 
assumed  phase  changes  from  SV  to  SH  to  P  etc.  This  occurs  because  the  propagation 
asimuth  for  SH  motion  is  set  to  be  +  or  -  90  *  from  the  particle  motion. 

For  the  second  example,  two  independent  samples  of  random  noise  were  combined 
on  the  N,  E  and  Z  components  of  motion  to  simulate  one  set  of  arrivals  traveling  at  0  * 
asimuth,  72  *  angle  of  incidence  with  an  amplitude  of  1000  and  another  set  of  arrivals 
traveling  N45*E,  55*  angle  of  incidence  and  amplitude  548  (amplitudes  are  relative). 
In  terms  of  the  independent  non-coherent  time  series  x(t)  and  y(t): 

N(t)=3.Qx(t)+y(t) 

E(t)=y(t)  (9) 

Z(t)=x(t)+y(t) 

Using  these  numbers  in  equation  (8),  we  obtain  values  of  6  *  for  the  asimuth  and 
67*  for  the  angle  of  incidence  (assuming  an  asimuth  of  6*  for  the  radial  trace).  The 
results  of  the  polarisation  analysis  are  shown  in  Figure  3.  The  asimuth  trace  is  centered 
about  6  *  over  the  length  of  the  trace  in  excellent  agreement  with  the  theory.  Error  for 
the  asimuth  is  large  and  is  due  to  the  fact  that  the  two  arrival  sets  are  traveling  at 
azimuths  45  *  apart;  the  worst  possible  case.  For  the  apparent  angle  of  incidence,  the 
measured  value  of  67  *  is  also  in  excellent  agreement  with  theory  and  in  this  case  the 
error  is  very  small.  The  small  error  arises  because  the  two  arrival  groups  are  traveling  at 
a  similar  angle  of  incidence  in  the  vertical  plane  at  6  *  asimuth. 

The  adaptive  polarization  analysis  program  is  a  useful  tool  for  locating  earthquakes 
when  only  a  single  three-component  data  set  is  available.  The  azimuth  from  the  initial 
P- wave,  where  presumably  the  error  will  be  smallest  and  no  conflicting  arrivals  are 
present,  combined  with  the  distance  obtained  from  the  S-P  time  gives  a  fairly  accurate 
location. 

Synthetic  Seismograms 

The  synthetic  seismograms  were  computed  using  the  locked  mode  method  (Harvey, 
1981)  at  frequencies  from  0  to  5  Hz.  Regional  geology  in  the  Northeast  required  that 
two  different  velocity/attenuation  models  be  used  for  the  computation;  a  New  England 
model  for  the  events  in  New  England  and  southern  New  York,  and  a  Grenville  model  for 
events  within  the  Grenville  province  (Figure  4,  Table  1).  The  New  England  velocity 
model  was  derived  by  Taylor  et  aJ.,  (1980)  using  regional  travel  times  of  P  and  S  waves 
recorded  across  the  Northeastern  United  States  Seismic  Network  operated  by  several 
academic,  governmental  and  private  institutions  in  the  northeast.  The  Grenville  model 
was  from  unpublished  refraction  results.  The  frequency  and  depth  dependent  Q  model 
for  the  New  England  model  was  adapted  from  work  by  Mitchell  (1981)  and  a  simple  fre¬ 
quency  dependent  Q  was  used  for  the  Grenville  model.  Although  there  are  significant 
differences  in  the  Q  models,  there  is  little  effect  on  the  synthetics  because  the  frequencies 
and  ranges  at  which  they  are  being  computed  are  relatively  small. 
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Figure  S.  Results  of  trying  to  trsck  two  different  sets  of 
random  noise  traveling  in  different  directions. 
One  set  is  traveling  to  the  northeast;  the  other  is 
traveling  due  north  and  is  almost  twice  the 
amplitude  of  the  first. 
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Figure  4.  Location  of  northeastern  events  and  stations 
used  in  this  study.  Events  are  indicated  by  cir¬ 
cles  and  a  date.  The  two  stations,  RSNY  and 
SRNY  are  indicated  by  triangles. 
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Synthetics  were  computed  for  each  event  at  several  depths  with  the  reported  focal 
mechanism  listed  in  Table  2.  The  instrument  response  was  convolved  with  the  synthet¬ 
ics  for  comparison  to  the  data. 

Station  RSNY  lies  within  the  Grenville  structure  close  to  the  boundary  with  the 
New  England  Structure.  While  the  travel  paths  to  RSNY  from  New  England  events 
must  pass  through  both  structures,  the  majority  of  the  path  is  within  the  New  England 
structure  and  so  the  New  England  structure  is  a  good  approximation  to  the  true  path. 

Throughout  this  report,  only  the  first  20  or  so  seconds  of  the  synthetic  depth  sec¬ 
tions  are  shown  as  we  wish  to  emphasise  the  characteristics  of  the  depth  phases  with 
focal  mechanism,  depth,  range,  and  velocity  model.  Some  of  the  phases  are  identified  in 
the  synthetic  sections  with  labeled  lines.  Reflections  off  of  the  13  km  discontinuity  in  the 
New  England  model  and  the  4  km  discontinuity  in  the  Grenville  model  are  indicated 
with  an  »  (e.g.  PiP).  Mantle  reflections  are  indicated  by  an  m.  The  depth  phases  are 
characterised  by  progressively  later  arrivals  with  increased  depth  and  the  primary  phases 
by  either  very  little  change  in  arrival  time  with  depth  for  the  direct  arrivals  or  progres¬ 
sively  earlier  arrivals  with  increased  depth.  In  some  of  the  sections,  the  depth  phases 
dominate,  while  in  others,  the  primary  phases  dominate.  This  is  a  function  of  the 
theoretical  radiation  pattern  predicted  by  the  input  focal  mechanism.  Because  the  input 
velocity  structure  will  influence  the  take-off  angle  of  each  phase  and  hence  its  position  on 
the  focal  sphere,  both  focal  mechanism  and  velocity  model  contribute  to  the  relative 
amplitudes  of  primary  and  depth  phases. 


TABLE  1 


New  England  Velocity/Q  Model  (Taylor  et  at. ,  1080) 


Thickness  (km) 

Vp  (km/sec) 

Vs  (km/sec) 

P  (gm/cm3) 

Q. 

9* 

2 

6.0 

1  ■ 

<xys 

MMjBVjga 

13 

6.1 

25 

7.0 

8.1 

Grenville  Velocity/Q  Model 

Thickness  (km) 

KMB9E51 

H3QI 

Q~ 

qjt 

4 

6.1 

3.5 

2.5 

1100+150  / 

5Qtf/9 

31 

6.6 

3.7 

2.7 

8.1 

4.6 

3.2 

The  success  of  determining  focal  depth  using  depth  phases  and  comparing  them  to 
synthetics  depends  largely  on: 

(1)  good  approximations  of  crustal  velocity  structure,  and 
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(2)  recovery  and  correct  identification  of  the  phases  in  the  data. 

Here  we  use  one  earthquake  to  illustrate  how  the  synthetic  seismograms  are  affected  by  2 
different  velocity  models  and  focal  mechanisms. 

The  earthquake  chosen  is  the  Gaza,  N.H.,  magnitude  4.5  event  of  January  19,  1982. 
Since  the  source  is  in  the  Appalachian  Province  and  the  receiver  (RSNY)  is  in  the  Gren¬ 
ville  Province  (see  Figure  4)  we  calculated  separate  synthetics  using  the  two  velocity 
models  (Table  1)  to  test  which  is  more  appropriate. 

Figure  5  shows  radial  and  vertical  components  for  the  initial  15  seconds  of  the  syn¬ 
thetic  seismograms  calculated  for  two  velocity  models.  RSNY  is  the  theoretical  receiver, 
a  distance  of  267  km  and  an  azimuth  of  295  *  from  the  New  Hampshire  source.  The 
velocity  models  and  Moho  reflections  are  shown  schematically  for  comparison.  The  focal 
mechanism  is  given  by:  strike  “  280  * ,  dip  *■  75  *  and  rake  ■»  -11  * .  The  input  focal 
depths  are  comparable. 

As  expected,  the  prominent  phases  and  their  arrival  times  are  different  for  the  two 
different  velocity  models.  The  Pn  phase  is  nodal  for  the  radiation  pattern,  despite  small 
differences  in  take-off  angles  for  the  two  models.  It  is  simply  indicated  at  the  appropri¬ 
ate  arrival  time.  For  the  New  England  model,  PmP  is  small  relative  to  tPmP.  They 
both  arrive  before  the  intra-crustal  reflections  PiP  and  aPiP.  For  depths  greater  than  10 
km,  however,  PiP  arrives  before  a  PmP.  PiP  and  aPiP  are  sharp  arrivals  because  of  the 
large  velocity  contrast  at  15  ’tn  in  the  crustal  model.  Vp  increases  from  6.1  km/sec  to 
7.0  km/sec  at  this  boundary 

For  the  Grenville  model,  Pg  precedes  the  Moho  reflections  (for  a  source  depth  of  7 
km).  The  amplitude  of  aPmP  relative  to  PmP  is  smaller  for  the  Grenville  model  than 
for  the  New  England  model  because  the  upgoing  aPmP  is  crossing  a  fairly  sharp  internal 
discontinuity  near  the  source  while  the  downgoing  PmP  only  crosses  this  boundary  once, 
near  the  receiver.  In  the  case  of  New  England  structure,  the  sharp  discontinuity  is  below 
the  source,  instead  of  above  it.  In  addition,  PmP  is  close  to  a  nodal  plane  and  small 
differences  in  take-off  angle,  resulting  from  differences  in  the  velocity  structures,  may  be 
contributing  to  the  difference  in  the  amplitudes  of  PmP. 

Even  though  the  two  velocity  models  produce  strikingly  different  synthetic  seismo¬ 
grams,  it  is  important  to  note  that  for  the  depths  and  models  shown,  the  time  between 
the  arrivals  of  the  downgoing  and  upgoing  Moho  reflections  is  the  same.  Thus,  if  one 
can  independently  identify  these  two  phases  in  the  actual  data,  a  match  (ignoring  abso¬ 
lute  travel  times)  to  the  Grenville  synthetic  would  suggest  a  depth  of  7  km  and  a  match 
to  the  New  England  synthetic  would  suggest  a  depth  of  6  km.  Given  this  time 
difference  between  the  two  arrivals,  the  choice  of  velocity  model  has  had  only  a  small 
effect  on  depth  determination.  If,  however,  one  uses  the  computed  phase  travel  times  to 
help  identify  depth  phases  in  the  real  data,  the  variation  in  travel  time  with  velocity 
model  could  lead  to  misinterpretations  of  the  real  data. 

Next,  in  Figure  6  we  show  15  seconds  of  synthetics  calculated  for  2  different  focal 
mechanisms.  The  mechanism  shown  at  the  bottom  was  determined  from  P- wave  first 
motion  data,  mainly  from  the  Northeast  United  States  Seismic  Network  (Pulli  et  ol . , 
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1083).  The  mechanism  at  the  top  of  the  figure  differs  mainly  in  the  quadrants  of  dilata-  9 

tion  and  compression.  The  asimuth  to  station  RSNY  is  rery  close  to  a  nodal  plane  in 
both  instances,  so  there  is  not  a  great  difference  in  radiation  pattern.  Thus,  the  same 
initial  phases  appear  on  both  sets  of  seismograms,  but  the  amplitudes  of  these  phases  are 
different.  If  theoretical  radiation  patterns  can  be  used  to  predict  actual  phases  in  the 
data,  then  it  will  be  helpful  to  have  independently  determined  fault  plane  solutions. 
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January  19,  1982,  a*  267  Km,  Azimuth*  295. 

Two  velocity  models  end  resulting  synthetic  seismograms. 

A)  New  England  Model,  source  d*pth—fl  km.  Source- to- receiver  Moho, 
reflections  shown  schematic  ally.  15  seconds  of  radial  and  vertical  com¬ 
ponent  synthetics,  with  phases  identified.  PiP  and  tPiP  are  internal 
reflections  at  the  15  km  boundary. 

B)  Grenville  model,  source  depth— 7  km.  Source* to- receiver  Moho, 
reflections  shown  schematically.  15  seconds  of  radial  and  vertical  com¬ 
ponent  synthetics  with  phases  identified. 
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January  19,  1982 

Gaza,  New  Hcrpshire 

Synthetics  for  2  Fault  Plane  Solutions 

New  England  Velocity  Model 

A*  267  Mn,  Azlmith=®5# 


Figure  0.  Two  fault  plane  solution*  and  resulting  synthetic 
seismograms.  Fifteen  seconds  of  the  radial  and 
vertical  components  are  shown  for  each  of  the 
two  different  input  solutions.  The  azimuth  to 
RSNY  is  indicated  by  a  line  on  the  fault  plane 
solution,  a  lower  hemisphere  projection  with 
compreasional  quadrants  filled  in. 
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SECTION  m 
DATA 


In  the  northeastern  United  States  and  Ontario,  Canada,  eight  local  and  regional 
earthquakes  and  one  quarry  blast  were  selected  for  study.  The  events  were  chosen  for 
their  size,  location,  and  availability  of  data.  A  list  of  these  events  is  given  in  Table  2 
along  with  location,  depth,  magnitude,  strike,  dip,  and  rake  of  the  focal  mechanism,  the 
region  (either  New  England  or  Grenville),  and  references  for  the  focal  parameters.  The 
event  and  station  locations  are  plotted  in  Figure  4,  showing  their  spatial  distribution. 
Many  of  the  earthquakes  in  Table  2  have  been  studied  by  others  and,  for  some,  the  focal 
depths  and  focal  mechanisms  are  well  determined  and  will  provide  an  excellent  base  for 
comparison. 

In  addition  to  reviewing  the  literature,  we  examined  phase  data  from  local  seismic 
networks  originally  used  to  invert  for  hypocentral  location.  In  three  cases  the  estimated 
error  on  focal  depth  was  so  large  that  to  use  the  reported  depth  would  be  misleading. 
For  these  we  indicate  (Table  2)  an  unknown  depth. 

Most  of  the  data  were  recorded  by  RSNY,  in  northern  New  York,  at  44.55 ‘N, 
74.53  *  W.  Several  of  the  recent  events  were  also  recorded  by  SRNY,  a  very-broad-band 
seismic  station  operated  by  Rondout  Associates  in  Stone  Ridge,  New  York,  at  41.85  *  N, 
74.15*  W. 

Data  for  the  Kuril/Kamchatka  portion  of  this  study  were  recorded  by  the  fourth 
Ocean  Sub-bottom  Seismometer  (OSS- IV)  emplaced  380  m  below  the  sea  floor  at  the 
bottom  of  DSDP  hole  581C  (44  *  N,  100  *  E).  The  data  spanned  the  period  from  Sep¬ 
tember  12  to  November  10,  1082  and  the  events  selected  for  study  are  listed  in  Table  3. 
The  closest  event  to  OSS-IV  in  574  km  away  and  the  furthest  1307  km.  The  large  dis¬ 
tances  make  the  computation  of  synthetic  seismograms  using  the  locked  mode  method 
unreasonable. 

Crosstalk  and  skew  problems  hamper  analyses  of  the  OSS-IV  data.  The  skew  was 
easily  corrected  but  the  amount  of  crosstalk  is  difficult  to  judge.  Theoretically  the 
crosstalk  problem  effects  only  the  azimuth  and  apparent  angle  of  incidence  determina¬ 
tions  and  not  the  isolation  of  the  depth  phases  needed  for  depth  determination.  Relative 
amplitudes  of  the  phases  may  not  be  correct,  but  the  polarization  state  filter  does  not 
preserve  absolute  amplitudes  anyway,  so  the  point  is  unimportant;  it  is  only  mentioned 
here  for  the  sake  of  completeness. 

A  map  of  the  event  locations  in  Table  3  and  OSS-IV  is  given  in  Figure  7.  The 
events  were  classified  as  being  either  excellent,  good  or  bad  for  depth  determination 
based  on  the  signal-to-noise  ratio  and  impulsiveness  of  the  records.  Recordings  that  were 
emergent  and  unlikely  to  show  clear  depth  phases  without  a  great  deal  of  processing  were 
classified  as  bad.  The  excellent  (*)  and  good  (x)  data  seem  to  be  associated  with  the 
downgoing  Pacific  plate  and  the  bad  (o)  data  with  the  shallow  events  within  the  accre¬ 
tionary  wedge  of  the  Kamchatka  peninsula.  The  number  of  events  studied  is  too  small 
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LIST  OF  EVENTS  AND  EVENT  PARAMETERS 
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events  used  in  this  study.  *  —  excellent  events, 
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to  attach  much  significance  to  this  observation;  however,  the  results  are  easily  explained 
by  the  tectonic  setting  of  the  region.  In  general,  energy  that  must  pass  through  the  tec¬ 
tonic  boundary  between  the  Asian  and  Pacific  plates  is  attenuated  much  more  than  that 
which  has  a  travel  path  wholly  in  the  Pacific  plate. 


TABLE  3 

Selected  Kuril/Kamchatka  Events  Recorded  By  OSS- IV 


Event 

Date 

Time 

Lat 

Lon 

Depth 

52.15 

47.10 


10/03/82 

10/11/82 

10/11/82 

10/11/82 

10/25/82 

10/25/82 

11/13/82 

11/15/82 

11/18/82 


21:20:44.0 

14:05:30.0 

23:10:25.5 

18:53:06.8 


00:32:57.2 

15:31:26.0 

15:33:22.5 

15:35:22.6 

12:05:47.3 

15:38:58.6 

10:58:46.0 

10:18:17.1 

00:51:12.6 


50.05 

54.84 

55.25 


154. 


158.68 

153.5 

158.10 

158.18 


161.37 

150.61 

150.55 

150.66 

157.2 

162.18 

157.60 

166.52 

163.76 


A  search  was  made  for  data  from  station  MAJO  in  Matsaahiro,  Japan  to  aid  in  the 
Kuril/Kamchatka  study.  Unfortunately,  only  long  period  data  was  available  for  the 
events  listed  in  Table  3. 
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SECTION  IV. 
DATA  ANALYSIS 


In  this  section  we  show  how  we  estimate  source  depths  of  the  earthquakes  presented 
in  the  previous  section.  First,  the  northeast  North  American  events  are  compared  to 
synthetic  seismogram  depth  sections  and  a  best-fit  depth  is  selected.  This  is  followed  by 
a  brief  presentation  of  the  Kuril/Kamchatka  OSS  data,  which  turns  out  to  be  ill-suited 
for  depth  determination. 

New  Hampshire,  January  10, 1982 

It  would  be  very  difficult  to  compare  the  raw  data  recorded  at  RSNY  to  the  syn¬ 
thetic  seismograms  (compare  the  top  of  Figure  8  with  the  synthetics  in  Figure  0).  With 
the  polarization  filtering,  however,  we  have  been  able  to  recover  polarised  phases  (bot¬ 
tom  of  Figure  8)  that  were  all  but  lost  in  the  P-coda  of  the  unfiltered  seismograms. 
After  polarization  filtering,  we  low-pass  filter  the  data  to  remove  frequencies  greater  than 
5  Hz,  which  are  not  modeled  by  the  synthetics. 

The  following  figures  show  comparisons  of  the  vertical-component  RSNY  processed 
data  to  the  vertical-component  synthetic  depth  section  calculated  for  the  New  England 
model.  Twenty  seconds  are  shown.  Ws -began  by  inserting  the  data  at  a  depth  of  about 
9  km  (Figure  10).  There  are  several  arrivals  in  the  data  that  match  predicted  arrivals. 
These  are,  in  order,  PmPt  »PmP ,  the  dual  arrivals  PiP  and  Pg,  and  finally  an  upgoing 
phase  arriving  toward  the  end  of  this  time  window.  Clearly,  the  match  is  not  perfect. 
Signal  arrives  in  the  data  later  than  the  appropriate  time  for  tPiP ,  and,  subsequently 
until  the  late  phase  (in  the  data,  about  14  seconds  after  the  initial  P),  there  is  not  a 
good  correlation  between  data  and  synthetics.  A  comparison  of  travel  times  for  the  first 
arriving  phase  in  the  data  and  synthetics  suggests  that  they  may  not  be  the  same  phase. 
The  travel  time  of  synthetic  PmP  at  a  source  depth  of  0  km  (with  the  given  input  velo¬ 
city  model)  is  40.40  seconds,  whereas  the  travel  time  of  the  first  phase  in  the  data  is 
38.50  seconds.  Allowing  for  uncertainties  in  the  actual  origin  time  of  the  earthquake  and 
for  uncertainties  in  crustal  structure,  they  could  be  the  same  phase.  If,  however,  we 
assume  that  the  origin  time  and  velocity  model  are  appropriate,  we  can  match  several 
early  P-phases  by  aligning  the  data  at  a  depth  of  approximately  4.5  km  as  illustrated  in 
Figure  11.  Since  the  real-data  travel  time  is  more  appropriate  for  Pn,  the  first  arrival  is 
now  aligned  with  the  predicted  Pn  arrival  time  (compare  Figures  10  and  11).  For  this 
shallower  depth,  the  data  conform  to  predicted  arrival  times  for  Pn,PmP,tPmP  and  for 
the  late  up-going  phase.  Again,  the  correlation  among  all  phases,  real  and  synthetic,  is 
not  perfect. 

At  this  point,  it  is  clear  that  a  great  deal  depends  on  the  correct  identification  of  the 
arrivals  in  the  data.  In  addition  to  travel  times,  the  azimuth  and  apparent  angles  of 
incidence  can  be  used  to  help  identify  phases.  In  Figure  12  we  show  the  results  of  an 
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Figure  8.  Seismograms  recorded  at  RSNY  of  the  January  10,  1082,  New  Hampshire 
earthquake,  A— 267  km.  The  top  three  traces  are  the  vertical,  north-south, 
and  east-west  components,  unfiltered.  The  bottom  three  traces  are  the  three 
components  after  state-filtering  for  rectilinear  particle  motion. 
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Figure  10.  Same  as  Figure  0  except  the  real  data  are 
included  in  the  depth  section  at  a  depth  of  about 
0  km.  See  text  for  disc  union. 
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Figure  11. 


Like  Figure  10  except  notice  that  the  reel  date 
are  plotted  at  a  depth  of  about  4.5  km.  See  text 
for  discussion. 
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Figure  It.  Adaptive  polarisation  of  state-filtered  seismo¬ 
gram  recorded  at  RSNY.  The  values  of  azimuth 
and  apparent  angle  of  incidence  are  given  with 
estimated  errors. 
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adaptive  polarization  analysis  of  the  first  20  seconds  of  the  same  state-filtered  seismo¬ 
grams  that  we  have  been  comparing  to  the  synthetics.  The  first  three  traces  '.a  the  figure 
are  the  vertical,  and  the  adaptive  radial  and  transverse  horizontal  components.  The 
direction  of  maximum  signal  strength  as  a  function  of  time  represented  by  the  azimuth 
and  apparent  angle  of  incidence,  is  indicated  (with  error  estimates)  below  the  RZ  pro¬ 
duct  trace  in  Figure  12.  The  positive  values  for  the  product  of  the  radial  and  vertical 
components  indicate  that  all  phases  shown  are  arriving  at  the  receiver  as  P-waves. 
Adaptive  polarization  analysis  of  the  synthetic  seismograms  (not  shown  here)  indicates 
the  same  thing;  thus,  at  the  very  least,  we  know  that  the  last  conversions  were  to  P- 
waves  for  both  the  real  and  synthetic  data. 

Interestingly,  we  noted  that  the  last  arrival,  which  had  matched  a  predicted  phase 
in  the  synthetics  at  both  the  4.5  and  9  km  depth  positions,  is  probably  not  even  modeled 
by  the  synthetics.  Though  highly  polarized  and  hence  a  prominent  phase  in  the  data,  it 
is  fourteen  degrees  off  the  correct  azimuth  from  the  source.  Unlike  the  two  other  out-of¬ 
line  azimuths  (see  Figure  12),  the  error  estimate  is  small  and  does  not  allow  enough 
uncertainty  to  put  it  in  line.  Furthermore,  this  late  phase  is  approaching  the  receiver  at 
a  shallower  angle  than  any  other  arrival  (see  the  angles  of  incidence),  yet  the  late  phase 
in  the  synthetic  seismograms  approaches  the  receiver  at  the  steepest  angle  in  the  same  20 
second  time  window.  A  steep,  late  arrival  would  be  consistent  with  a  multiply-reflected 
phase. 


In  the  data,  the  phase  immediately  preceeding  the  last  phase  (by  almost  1  second)  is 
more  consistent  with  the  synthetics  because  it  approaches  at  a  steeper  angle.  Using  that 
information,  a  different  match  can  be  obtained  by  inserting  the  data  at  a  depth  of  6  km 
(Figure  13).  Here,  the  first  arrival  is  a  predicted  PmP  and  there  are  phases  matching 
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the  predicted  sPmP  and  the  dual  arrival  of  Pg  and  PiP,  as  well  a a  the  late  multiply- 
reflected  phase.  This  interpretation,  however,  gives  the  largest  mismatch  in  travel  times. 

One  can  see  that  it  is  not  easy  to  choose  among  the  above  interpretations  ,  because 
several  are  possible  yet  none  provides  a  perfect  match.  The  difficulties  of  identifying  and 
correlating  phases  from  a  single  set  of  three-component  data  are  well-illustrated  by  this 
example.  Evidently,  this  particular  earthquake  is  difficult  to  pin  down  as  the  range  of 
estimates  for  its  focal  depth  is  3-11  km,  using  a  variety  of  approaches  (see  Table  4). 


Adirondack  Mountains,  August  31, 1033 

No  depth  has  been  reported  in  the  literature  for  the  31  August,  1082,  Adirondack 
event.  Its  location  is  153  km  from  RSNY  at  a  station  to  event  asimuth  of  170  * .  Depth 
phases  would  be  very  difficult  to  identify  if  just  the  band  pass  filtered  data  rotated  to  the 
vertical,  radial,  and  transverse  directions  of  motion  were  examined  (Figure  14,  top).  In 
this  case,  the  polarization  state  Alter  was  used  effectively  to  pass  only  rectilinear  motion 
and  help  in  the  depth  phase  identification  (Figure  14,  bottom). 

The  synthetic  seismogram  depth  section  computed  for  this  event  exhibits  a  strong 
tPmP  phase,  well  defined  SmP,  pSmP,  and  tPmS  phases  and  a  weak  Pg  phase  (Figure 

15) .  The  strike,  dip,  and  rake  for  this  event  have  not  been  published  and  a  generic  fault 
plane  solution  for  this  area  based  on  other  Adirondack  events  was  used  in  the  synthetic 
program.  Thus,  the  relative  amplitudes  of  the  phases  do  not  match  the  data  (Figure 

16) .  The  state  filtered  data  trace  in  Figure  16  has  been  aligned  with  the  synthetics  in 
absolute  time.  In  other  words,  no  movement  of  the  data  trace  with  respect  to  the  syn¬ 
thetics  along  the  time  axis  was  allowed.  With  this  constraint,  the  large  first  arrival  in 
the  data  trace  is  identified  as  the  Pg  arrival  and  the  second  arrival  is  the  $PmP  arrival. 
The  small  arrival  at  about  10  seconds  on  the  data  trace  may  be  the  pSmP  arrival. 

The  beat  fit  of  the  data  to  the  synthetics  is  at  7  km  depth.  We  note  however,  that 
if  amplitudes  are  ignored  (as  they  should  be  in  this  case)  the  fit  of  the  data  to  the  syn¬ 
thetics  is  nearly  as  good  for  any  depth  between  7  and  13  km. 


Maine,  May  20,  1083 

The  reported  hypocenter  of  the  Maine  event  is  at  44.50  *  N,  70.41  *  W,  at  a  depth  of 
2  km.  This  location  is  328  km  due  east  of  RSNY,  the  largest  distance  to  any  of  the 
northeast  U.S.  events  studied.  The  data  are  characterized  by  a  small  Pn  phase  arriving 
46  seconds  after  the  rupture  followed  5  seconds  later  by  a  rather  large  emergent  Pg.  The 
coda  amplitude  remains  fairly  constant  from  the  Pg  arrival  to  the  Lg  arrival  38  seconds 
later  (Figure  17).  Synthetic  seismograms  computed  for  this  event  show  large  amplitude 
•PmP,  PiP  and  »PiP  phases  and  several  unidentified  depth  phases  (Figure  18).  The 
coda  of  the  Pg  arrival  in  the  data  is  not  well  matched  by  the  synthetics  for  a  2  km 
source  depth.  Of  more  importance  is  the  mismatch  of  amplitudes  of  the  tPmP  arrival. 
The  tPmP  arrival  in  the  data  has  the  same  small  amplitude  as  the  PmP  arrival.  A 
better  fit  of  the  data  to  the  synthetics  occurs  at  12  km  source  depth  (Figure  19).  At  this 
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Figure  14.  Vertical,  radial,  and  transverse  components  of 
motion  for  the  August  31,  1982,  Adirondack 
earthquake  as  recorded  at  station  RSNY,  A— 156 
km.  Lower  traces  have  been  polarisation  state 
filtered  for  rectilinear  motion  to  enhance  the 
depth  phases. 


Figure  10.  Best  fit  of  the  RSNY  data  to  the  synthetics  for 
the  Adirondack  earthquake  (7  km).  There  has 
been  no  published  depth  for  this  event.  See  Fig¬ 
ure  15  for  details  of  the  synthetics. 
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Figure  18.  Synthetic  seismogram  depth  section  for  the 
Maine  earthquake  A— 328  km,  strike— 355*, 
dip— 52  * ,  rake— 78  * ,  station  to  event 

aiimuth— 90  * . 
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depth  the  aPmP  phase  arrives  in  the  Pg  coda  where  the  amplitudes  are  better  matched 
and  the  arrival  times  of  the  distinct  phases  in  the  data  are  modeled  well  by  the  synthet¬ 
ics.  Regional  network  data  from  36  stations  placed  the  depth  of  this  event  at  1.8  km. 
Support  is  lent  to  this  depth  by  an  aftershock  survey  conducted  shortly  after  the  main 
shock  that  obtained  a  well  determined  depth  of  2.4  km  for  one  of  the  aftershocks  (Ebel 
and  McCaffrey,  1084). 


Goodnow,  New  York,  October  7, 1083 

The  Goodnow  sequence  of  earthquakes  occurred  in  the  central  Adirondack  moun¬ 
tains  of  New  York  on  October  7,  1983.  Following  the  mt  5.3  main  shock,  there  were 
nearly  100  aftershocks  with  magnitude  3.5  and  lower  that  occurred  over  the  next  two 
months.  Station  RSNY,  just  60  km  from  the  event  recorded  the  P-waves  reliably  but 
clipped  in  the  Lg  portion  of  the  wavetrain.  The  strike,  dip,  and  rake  of  the  main  shock 
from  first  motion  data  are  173  * ,  60* ,  and  110* .  Station  RSNY  is  at  an  azimuth  of 
162  *  and  lies  close  to  the  intersection  of  the  two  nodal  planes  of  the  main  shock.  Thus, 
the  record  character  for  synthetics  at  RSNY  is  sensitive  to  small  changes  in  focal 
mechanism. 

The  data  shown  in  Figure  20  have  been  rotated  to  the  vertical,  radial,  and 
transverse  directions  and  has  been  band-pass  filtered  between  .5  and  5  Hz  for  comparison 
to  the  synthetics.  Several  distinct  phases  are  evident  in  the  data. 

Figure  21  shows  the  synthetic  seismogram  depth  section  for  this  event.  The  Pg 
arrival,  strong  when  the  source  is  deep,  weakens  significantly  with  foci  near  to  the  sur¬ 
face.  The  predominant  phases  are  the  «PtP  and  the  SmP  which  exhibit  strong  ampli¬ 
tudes  throughout  the  section. 

A  match  of  the  data  to  the  synthetics  is  shown  in  Figure  22.  The  initial  Pg  arrival 
in  the  data  trace  is  about  .5  seconds  earlier  than  the  same  phase  in  the  synthetic  section 
and  the  second  arrival  is  not  matched  at  all  by  the  synthetics.  Beyond  that  however,  the 
tPcP,  aScP,  pPmP  and  SmP  phases  are  all  matched  well  for  a  source  depth  of  10  km. 
Despite  the  sensitivity  of  the  synthetics  to  small  changes  in  the  focal  mechanism  at  this 
azimuth,  the  relative  amplitudes  of  the  phases  are  fairly  well  matched  also. 

Two  depths  have  been  reported  for  this  event,  one  at  7  km  ±  1  and  one  at  14  km  ± 
5.  The  shallow  depth  is  based  on  aftershocks  and  on  teleseismic  observations  of  the  main 
event,  while  the  deeper  depth  is  based  on  local  travel  times.  The  depth  determined  here 
lies  between  these  two  estimates,  range  of  depths. 

The  large  size  of  the  Goodnow  event  (mt  =  5.3)  raises  the  possibility  of  a  complex 
source  function.  If  this  were  the  case,  the  complexities  could  be  falsely  identified  as 
depth  phases  and  an  improper  depth  determined.  At  this  point,  we  do  not  know 
whether  or  not  source  complexities  have  influenced  the  depth  for  the  Goodnow  event. 
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Figure  21.  Synthetic  seismogram  depth  section  for  the 
Goodnow  earthquake.  A  —69  km,  strike— 173  * , 
dip— 60',  rake— 110*,  station  to  event 

azimuth— 163  * . 
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Figure  22.  Best  fit  of  the  RSNY  data  to  the  synthetics  for 
the  Goodnow  earthquake  (10  km).  The  pub¬ 
lished  depths  range  between  7  and  12  km.  See 
Figure  21  for  details  of  the  synthetics. 
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Ontario,  October  11, 1983 

The  Ontario  event  occurred  123  km  from  station  RSNY.  It  was  a  magnitude  4.2 
event  and  had  a  reported  depth  of  15  km.  Although  the  travel  path  between  the  event 
and  RSNY  was  wholly  within  the  Grenville  province  (Figure  4),  distinct  phase  are 
difficult  to  identify  in  the  data.  Figure  23  shows  the  data  rotated  to  the  vertical,  radial, 
and  transverse  directions.  From  the  initial  Pg  arrival,  the  coda  gradually  increases  in 
amplitude  for  the  next  4  seconds. 

Synthetic  seismograms  for  this  event  (Figure  24)  have  a  small  Pg  phase  followed  by 
a  large  amplitude  PmP  phases  a  few  seconds  later.  Well  into  the  coda  there  is  a  well 
defined  set  of  arrivals  consisting  of  the  SmP,  PmS ,  pSmP,  and  tPmS  phases 

The  fit  of  the  data  to  the  synthetics  is  poor  (Figure  25).  Neither  the  arrival  times 
nor  the  amplitude  characteristics  are  matched.  The  depth  that  we  have  assigned  to  this 
event  (14  km)  is  based  primarily  upon  the  spacing  of  the  (presumed)  Pg  and  PmP 
phases. 


Adirondack,  New  York,  October  23, 1984 

On  October  23,  1984,  a  magnitude  3.4  earthquake  occurred  in  the  southwestern 
Adirondacks,  near  Johnsburg,  New  York.  Because  it  was  recorded  by  RSNY  and  SRNY, 
it  afforded  us  the  chance  to  model  seismograms  at  two  different  receivers  and  thus 
increase  our  confidence  in  the  interpretations.  We  estimate  a  source  depth  of  12  km  and 
a  strike,  dip,  and  rake  of  353  * ,  60  * ,  and  17  *  for  this  earthquake  based  on  the  results 
described  below.  Neither  the  depth  nor  the  fault  plane  solution  could  be  determined  by 
local  network  data,  though  many  stations  recorded  the  event  and  its  epicentral  location 
is  fairly  well  constrained  with  an  azimuthal  gap  of  43  * .  The  nearest  network  station, 
however,  is  52  km  from  the  source. 

Synthetics  to  both  RSNY  and  SRNY  were  calculated  with  the  Adirondack  velocity 
model.  Initially,  we  tried  the  fault-plane  solution  of  the  Goodnow  earthquake  because  it 
is  well  constrained  and  is  similar  to  many  other  Adirondack  solutions  (e.g.  Yang  and 
Aggarwal,  1981).  Resulting  vertical  and  radial  synthetic  depth  sections  are  shown  in 
Figures  26  and  27  along  with  RSNY  band-pass  filtered  data  plotted  at  a  depth  of  12  km. 
The  data  and  synthetics,  though  aligned  at  the  same  time  with  respect  to  the  origin 
time,  do  not  look  much  alike.  The  lack  of  fit  between  data  and  synthetics,  particularly 
for  the  relative  amplitudes  of  different  phases,  suggests  that  the  fault  plane  solution  and 
hence  radiation  pattern  of  the  synthetics  is  incorrect.  A  much  better  match  is  achieved 
using  a  different  fault  plane  solution  of  strike  —  353  * ,  dip  -  60 ' ,  rake  —  17  *  (Figures 
28  through  31).  This  mechanism  has  a  greater  strike-slip  component  and  is  not  unlike 
several  previously  determined  Adirondack  mechanisms  (Schlesinger  -  Miller,  et  al., 
1983).  In  Figures  28  and  29  RSNY  vertical  and  radial  seismograms  are  again  plotted  at 
a  depth  of  12  km  where  the  arrival  times  and  relative  amplitudes,  especially  of  PmP  and 
•PmP ,  are  well  modeled  by  the  synthetics.  Furthermore,  the  polarity  of  the  3-cycle 
waveform,  tPmP,  is  opposite  that  of  PmP  on  both  data  and  synthetics.  Pg  is  nearly 
nodal  for  the  theoretical  seismograms  and,  in  fact,  no  signal  above  noise  was  detected  on 
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Figure  Bj.  Synthetic  seismogram  depth  section  for  the 
Ontario  earthquake.  A ->123  km,  strike— 71  * , 


seconds 


Figure  t5.  Best  fit  of  the  RSNY  data  to  the  synthetics  for 
the  Ontario  earthquake  (13  km).  The  published 
depth  for  this  event  is  15  km.  See  Figure  24  for 
details  of  the  synthetics. 
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figure  SO.  Vertical-component  synthetic  seismogram  depth  section  for  the  October  23, 
1984,  earthquake  with  the  RSNY  vertical  trace  plotted  at  12  km  depth. 
Synthetics  are  low-pass  Altered  from  3  Hz.  Data  are  band-pass  Altered,  1-3 
Hs.  Strike-173’,  dip-80’,  rake-110’,  A-115  km,  azimuth-336*. 
Receiver  response— RSNY  short  period. 
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Figure  27.  Radial-component  synthetic  seismogram  depth  section  for  the  October  23, 
1084,  earthquake  with  the  RSNY  radial  trace  plotted  at  12  km  depth.  Syn¬ 
thetics  are  low-pass  filtered  from  3  Hz.  Data  are  band-pass  filtered,  1-3  Hz. 
Strike— 173  * ,  dip— 80  * ,  rake— 110  * ,  A— 115  km,  azimuth- 336  * .  Receiver 
response— RSNY  short  period. 
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Figure  £8.  Vertical-component  synthetic  seismogram  depth  section  for  the  October  23, 
1084,  earthquake  with  the  RSNY  vertical  trace  plotted  at  12  km  depth. 
Synthetics  are  low-pass  filtered  from  3  Hz.  Data  are  band-pass  filtered,  1-3 
Hz.  Strike-353*,  dip-60',  rake-17*,  A-115  km,  azimuth-336*. 
Receiver  response— RSNY  short  period. 
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Figure  SO.  Radial-component  synthetic  seismogram  depth  section  for  the  October  23, 
1984,  earthquake  with  the  RSNY  radial  trace  plotted  at  12  km  depth.  Syn¬ 
thetics  are  low-pass  filtered  from  3  Hs.  Data  are  band-pass  filtered,  1-3  Hz. 
Strike— 353*,  dip— 80*,  rake— 17*,  A— 115  km,  asimuth— 338  * .  Receiver 
response— RSNY  short  period. 
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Figure  SO.  Vertical-component  synthetic  seismogram  depth  section  for  the  October  23, 
1984,  earthquake  with  the  SRNY  vertical  trace  plotted  at  12  km  depth. 
Synthetics  are  low-pass  filtered  from  3  Hs.  Data  are  band-pass  filtered,  1-3 
Hs.  Strike-353*,  dip-60*,  rake-17*,  A-195  km,  azimuth- 185  * . 
Receiver  response— SRNY  broad  band. 


Figure  SI.  Radial-component  synthetic  seismogram  depth  section  for  the  October  23, 
1984,  earthquake  with  the  SRNY  radial  trace  plotted  at  12  km  depth.  Syn¬ 
thetics  are  low-pass  filtered  from  3  Hz.  Data  are  band-pass  filtered,  1-3  Hz. 
Strike-353*,  dip-80*,  rake-17',  A-195  km,  azimuth- 1 85 * .  Receiver 
response— SRNY  broad  band. 
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the  observed  seismogram.  We  are  fairly  confident  of  our  interpretation  of  the  first  O 

arrival  at  RSNY  as  PmP  because  its  travel  time  and  phase  velocity,  6.8  km/sec,  are  both 
appropriate  for  PmP.  Phase  velocities  are  estimated  using  a  Vp  of  5.75  km/sec  at  the 
receiver  (as  reported  in  Taylor  and  Qualheim,  1983)  and  measuring  observed  angles  of 
incidence. 

Synthetics  for  the  second  fault  plane  solution  also  compare  reasonably  well  with  the  f? 

data  at  SRNY  for  the  same  depth  of  12  km  (shown  in  Figures  30  and  31).  Here  the 
travel  times  and  relative  amplitudes  of  Pg  and  sPmP  are  closely  matched.  On  the  other 
hand,  PmP  is  not  pronounced  in  the  data.  We  do  not,  however,  expect  the  fit  to  be  as 
good  because  SRNY  is  farther  from  the  source  and  lies  in  a  different  geologic  province. 
Nonetheless,  the  results  for  both  stations  support  the  interpretation  of  the  observed  $ 

phases  and  a  depth  estimate  of  12  km. 


Mid-Hudson  New  York  Earthquake  and  Quarry  Blast 

On  October  30,  1985,  a  small  earthquake  occurred  near  Amsterdam,  N.Y.,  only  4 
km  surface  distance  from  an  active  quarry.  This  gave  us  a  chance  (rare  in  the  eastern 
U.S.)  to  compare  an  earthquake  to  a  blast  source  from  nearly  the  same  location.  We 
attempted  to  model  the  earthquake  and  a  blast  to  station  SRNY,  a  distance  of  120  km 
and  117  km  from  the  respective  sources.  Several  lines  of  evidence  suggest  the  earthquake 
was  shallow,  but  we  caution  against  complete  confidence  in  a  one-station  interpretation 
because: 

1.  the  P- wave  signal-to-noiae  ratio  at  SRNY  is  poor  and 

2.  both  the  focal  mechanism  and  the  velocity  model  were  unknowns. 

A  visual  comparison  of  three-component  seismograms  (0.5-30  Hz)  for  the  earthquake 
and  for  a  large  quarry  blast  that  occurred  earlier  in  October,  1985  (Figure  32),  shows  a 
difference  in  P  and  S/Lg  coda  shape,  even  though  the  records  are  obviously  noisy.  One- 
to-five  Hz  band-pass  filtered  data  are  shown  in  Figure  33  to  illustrate  that,  despite 
differences  in  the  seismograms,  there  appear  to  be  5  and  Lg  phases  for  both  events. 
Unfortunately,  we  don’t  know  the  origin  time  of  this  particular  blast,  and,  because  the 
first  JP-wave  arrivals  are  so  unclear  for  both  events,  it  is  not  obvious  how  the  blast 
should  be  aligned  relative  to  the  earthquake  or  to  synthetics.  The  time  alignment  illus¬ 
trated  in  the  Figure  is  based  on  a  visual  match  of  large  amplitudes  on  the  vertical  com¬ 
ponents,  and  it  looks  sensible,  but  is  it  meaningful? 

We  subsequently  examined  digital  seismograms  of  network  stations  run  by 
Woodward-Clyde  Consultants.  Because  those  records  also  showed  similarities  in  relative 
timing  of  arrivals  between  the  earthquake  and  blasts  at  the  same  receiver,  we  concluded 
that  blast  phases  and  earthquake  phases  have  the  same  travel-times.  We  still  do  not 
understand  what  an  "5-wave"  from  the  blast  actually  is,  but  with  the  3-component 
SRNY  data  we  were  able  to  observe  the  particle  motion  and  neither  the  earthquake  nor 
the  blast  exhibit  5-type  motion  for  the  relatively  strong  arrival  that  has  a  group  velocity 
of  3.6  km/sec.  Nevertheless,  the  interpretation  of  "matching  phases"  shown  by  the 
alignment  of  the  time  series  in  Figure  33  is  what  we  will  use  to  compare  SRNY  data  to 
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Figure  85.  Verticil-component  synthetic  depth  section  compared  to  state-filtered 
seismograms  recorded  at  SRNY.  Top  trace  is  quarry  blast.  Third  from  top 
is  the  earthquake.  For  details  of  the  synthetics,  see  Figure  34. 
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synthetics.  Notice  that  the  first  visible  P-w&ve  is  later  on  the  earthquake  and  has  a 
group  velocity  of  5.6  km/sec.  The  earthquake’s  first  P-pick  at  a  nearby  regional  station 
(WPNY  of  Lamont-Doherty  Geological  Observatory’s  network),  however,  yields  a  group 
velocity  of  6.1  km/sec,  suggesting  that  an  earlier  arrival  at  SRNY  may  be  masked  by 
noise.  Polarization  filtering  of  SRNY  data  (as  seen  in  Figure  35  with  synthetics) 
improves  the  strength  of  rectilinearly  polarized  signal,  but  doesn’t  reveal  a  clear  P- 
arrival  earlier  than  the  one  with  a  group  velocity  of  5.6  km/sec. 

A  depth-section  of  vertical  synthetic  seismograms  is  shown  in  Figure  34.  The  sur¬ 
face  (actually  0.1  km)  depth  was  modeled  as  an  explosion,  while  the  remaining  depths 
were  modeled  as  a  double-couple  mechanism  with  strike  =  018  * ,  dip  =  60  ’ ,  rake  — 
156  * .  As  mentioned  above,  no  fault-plane  solution  could  be  determined  from  network 
first  motions  for  the  30  October,  1985  event,  so  we  used  a  generic  mid-Hudson  solution 
that  is  based  on  well-constrained  data  for  earthquakes  approximately  34  km  south  of 
Amsterdam,  N.Y.  (Houlday  et  al. ,  1984).  Also,  we  don’t  have  a  velocity  model  appropri¬ 
ate  for  travel  paths  between  Amsterdam  and  Stone  Ridge,  N.Y.  We  used  the  New  Eng¬ 
land  model,  rather  than  the  Adirondack  model,  because  geologically  it  is  an  Appalachian 
path. 

Prominent  early  arrivals  on  the  synthetic  seismograms  are  Pg,  PiP,  aPiP ,  and 
aPmP.  Though  the  synthetics  do  not  match  the  data  in  detail,  a  plausible  interpretation 
is  presented  in  Figure  35.  The  blast,  which,  we  verified  with  the  quarry,  must  be  a  stir- 
face  source,  and  the  travel  times  we  interpreted  in  Figure  33  are  close  to  those  of  the  ini¬ 
tial  two  P-pulses  modeled  as  a  surface  explosion.  As  expected,  the  blast  phase  with 
estimated  group  velocity  *  3.6  km/sec  is  not  accounted  for  by  theory.  The  first  detect¬ 
able  P-wave  from  the  earthquake  can  be  interpreted  as  a  depth  phase,  given  its  group 
velocity.  If  it  is  aPmP ,  then  the  earthquake  is  very  shallow.  If  it  is  aPiP ,  then,  from  its 
travel  times,  the  quake  is  about  7  km  deep.  We  prefer  the  shallow  depth  interpretation 
(shown  in  Figure  35)  for  at  least  two  reasons.  First,  the  nearest  network  station,  ROTD, 
has  a  2  Hz,  large-amplitude  arrival  with  a  group  velocity  of  approximately  2.8  km/sec. 
It  is  probably  a  fundamental-mode  Rayleigh  wave,  indicative  of  a  shallow  source. 
Second,  the  depth  estimate  given  by  the  network  hypocentral  location  is  4  ±  44  km,  yet 
the  nearest  station  is  only  20  km  away.  Such  runaway  error  estimates  can  arise  when 
the  true  depth  is  small  relative  to  epicentral  distance  to  the  nearest  station.  Finally,  the 
earthquake  (Me  —  2.7)  was  felt  locally  by  many  people  (Nottis,  personal  communica¬ 
tion),  suggesting  a  shallow  source.  Admittedly,  more  details  of  intensity  effects  would  be 
required  to  discriminate  a  very  shallow  foc*:a  from  a  7  km  focus. 

There  is  enough  evidence  to  conclude  that  the  Amsterdam  earthquake  is  shallow, 
but  given  uncertainties  in  faulting  mechanism  and  velocity  model,  along  with  a  poor 
signal-to-noise  ratio  at  SRNY,  the  identification  of  a  depth  phase  is  a  bit  shaky. 

Ardsley,  New  York,  October  19,  1985 

The  Ardsley  earthquake  (M  «  4.0)  of  October  19,  1985  occurred  in  southern  New 
York  state  at  40.98  *  N,  73.83  *  W  at  a  depth  of  5  km  (Figure  4).  It  was  felt  in  many 
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areas  of  New  York,  New  Jersey,  Connecticut,  as  well  aa  Pennsylvania,  Massachusetts, 
and  Rhode  Island.  The  largest  aftershock  (A i  **  3.0)  occurred  two  days  later,  also  at  a 
depth  of  5  km.  Both  the  main  shock  and  the  aftershock  were  recorded  by  station  SRNY, 
100  km  away  at  41.85  *  N,  74.15  *  W.  SRNY  is  a  three-component  broadband  digital  sta¬ 
tion  with  a  sample  rate  of  60  samples/second. 

A  comparison  of  the  main  shock  and  aftershock  seismograms  from  station  SRNY 
(Figure  36)  shows  a  remarkable  similarity.  The  traces  have  been  lined  up  with  the  first 
motion  of  the  aftershock  occurring  .72  seconds  after  the  first  motion  of  the  main  shock. 
The  traces  match  nearly  wiggle  for  wiggle  on  all  three  components.  The  aftershock 
matches  so  well  with  arrivals  well  into  the  main  shock  coda  that  the  main  shock  must  be 
a  complex  event  with  at  least  one  large  impulse  occurring  after  the  initial  impulse. 

It  is  well  known  that  the  rupture  history  of  a  large  seismic  event  can  be  extremely 
complex.  Studies  of  teleseismie  body  waves  have  identified  multiple  ruptures  (e.g. 
Kanamori  and  Stewart,  1078;  Trifunae  and  Brune,  1070),  Doppler  shifts  (Douglas  et  al. , 
1081),  body  wave  directivity  functions  (Khattri,  1078)  and  other  phenomena  as  contribu¬ 
tors  to  the  complexity  of  the  wavetrain.  In  contrast,  there  are  relatively  few  studies 
(Bakun  et  al. ,  1078)  that  have  identified  these  same  complexities  in  smaller  events.  In 
fact,  much  research  has  been  devoted  to  understanding  the  complexity  of  large  events  by 
assuming  that  small  events  in  the  same  area,  usually  aftershocks  of  the  large  event,  were 
generated  by  simple  point  sources  (e.g.  Hart  sell,  1078;  Kanamori,  1070;  Frankel  and 
Kanamori,  1083;  Mueller,  1085). 

We  have  capitalized  on  the  point-like  source  behavior  below  12  Hz  of  the  aftershock 
from  the  10  October  1085  Ardsley  earthquake  to  deconvolve  the  main  shock  using  the 
aftershock  as  an  empirical  Green’s  function.  By  asstiming  that  the  aftershock  is  a  point- 
source  Green’s  function,  the  effect  of  propagation  path,  recording-site  response,  and 
instrument  response  are  eliminated  as  sources  of  nutjor  error. 

The  method  of  using  small  aftershocks  as  empirical  Green’s  functions  to  deconvolve 
larger  events  is  based  on  the  assumption  that  the  smaller  event  is  a  simple  source  com¬ 
pared  to  the  larger  event.  The  closer  the  small  event  source  function  is  to  an  impulse  in 
time  and  space,  the  closer  the  event  approximates  the  true  Green’s  function. 

The  main  shock  displacement  can  be  written  as  a  convolution  of  the  source,  path, 
receiver,  and  instrument  responses  at  a  particular  source/receiver  geometry. 

D(t)  -  S(t)  *  P(t)  *  R(t)  *  I(t)  (10) 

where  *  denotes  convolution.  The  Green’s  function  for  the  same  source /receiver 
geometry  is  written: 

G(t)  =  P(t)  •  R(t)  *  1(1)  (11) 

Thus,  by  deconvolving  (10)  with  (11),  S(t),  the  far  field  displacement  time  history  for 
the  source  in  (10)  is  obtained. 
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The  information  gained  from  a  single  station  can  be  used  in  conjunction  with  the 
focal  fn**heniam  to  place  some  broad  constraints  upon  the  fault  dimensions.  To  obtain 
a  complete  record  of  the  large  fault  displacement  time  history,  data  from  at  least  three 
stations  are  needed.  From  the  time  difference  between  the  first  and  last  asperity  arrivals, 
At^  a  maximum  fault  dimension  (L _ )  can  be  determined  by  assuming  fault  propaga¬ 

tion  in  the  direction  of  the  receiver.  The  minimum  fault  dimension  (Lmln)  is  found  from 
the  maximum  time  difference  between  two  successive  asperity  arrivals  At _ 


L 


min 


(12) 


if  failure  propagation  is  less  than  5- wave  velocity  the  use  of  equations  (12)  will  give 
larger  than  the  actual  dimension. 

The  assumptions  made  in  this  analysis  purposely  oversimplify  the  faulting  mechan¬ 
ism  so  that  a  simple  solution  may  be  obtained.  The  results  presented  in  this  report  using 
the  above  methods  agree  well  with  the  theoretical  fault  dimensions  calculated  using 
empirical  formulas  for  mt  vs  fault  dimension  derived  by  Kanamori  and  Anderson  (1075), 
and  with  the  fault  dimension  determined  by  aftershock  studies. 

Very  little  processing  was  necessary  to  obtain  excellent  results  for  the  Ardsley 
events.  The  data  were  band  pass  filtered  between  0.5  and  12  Hz  with  a  phase-  free  filter 
and  then  the  horizontal  components  were  rotated  to  the  radial  and  transverse  directions 
for  station  SRNY  before  deconvolving.  The  data  were  also  processed  using  the  adaptive 
polarisation  analysis  method  and  the  polarization  state  filter  and  gave  valuable  informa¬ 
tion  about  the  polarisation  stability  of  the  coda.  However,  the  additional  polarization 
processing  did  not  change  the  results. 

In  Figure  37  we  show  the  result  of  deconvolving  the  first  few  seconds  of  P- wave  of 
the  Ardsley  earthquake  using  the  21  October  aftershock.  The  top  trace  in  the  plot  is  the 
signal  from  the  main  shock,  the  middle  trace  is  from  the  aftershock  and  the  bottom  trace 
is  the  deconvolution  result.  The  alignment  of  the  main  shock  and  aftershock  traces  in 
the  plots  is  shown  as  a  marker  in  the  deconvolution  result. 

There  are  three  distinct  peaks  in  this  result  separated  by  times  of  .57  and  .15 
seconds.  The  first  peak  (at  -0.75  seconds)  corresponds  to  the  alignment  of  the  first 
arrivals  from  both  events  and  the  third  peak  corresponds  to  the  best  visual  alignment  of 
the  traces  shown  in  Figure  36.  The  second  peak  cannot  be  easily  seen  by  aligning  the 
traces  visually. 

Because  only  one  station  was  used  in  the  analysis  presented  here,  relative  locations 
of  the  three  asperities  cannot  be  determined.  However,  using  equations  (12)  we  obtain 
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Figure  SI.  Deconvolution  of  the  Ardaley  m^in  shock  by  its 
largest  aftershock.  The  first  two  arrivals  (Pg 
&nd  PiP)  as  well  as  the  depth  phase  tPiP  are 
included  in  the  deconvolution  operator. 
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values  of  1.3  and  5.0  km  for  the  minifnmw  and  mMrimmw  fault  dimensions  namiming 
fault  propagation  at  the  shear  wave  velocity.  Better  estimates  of  the  fault  dimension 
limits  can  be  obtained  by  using  the  focal  mechanism  and  that  the  asperities  lie 

in  a  line  along  the  direction  of  the  rake.  The  strike,  dip,  and  rake  of  the  event  were 
reported  as  22  * ,  00  * ,  and  180  *  respectively  by  the  Lamont  Doherty  Geological  Obser¬ 
vatory  (Seeber,  personal  communication).  Station  SRNY  is  0.0  degrees  away  at  an 
azimuth  of  -15.4  * .  At  this  distance,  the  wavefront  plane  dips  at  about  60  * .  In  Figure 
38,  the  geometry  of  the  fault  and  wavefront  planes  are  shown.  We  have  assumed  that 
the  asperities  must  all  lie  along  the  line  of  the  rake  and  have  let  the  wavefront  plane  con¬ 
tain  the  first  asperity.  Thus,  the  distance  from  the  first  asperity  to  the  others  and  then 
to  the  wavefront  plane  divided  by  the  appropriate  velocities  must  equal  the  time 
differences  between  spikes  on  the  deconvolution  result.  From  Figure  38  and  equations 
(12)  the  minimum  and  maximum  fault  dimensions  are  given  by: 


^min/max 


_ ^num/to* _ 

V1  ±  coe(37.4°)  cos(30°)  Vp~x 


(13) 


Using  Vp  -  6  km/sec  ,  Vs  =  3.46  km/sec,  and  At_..  —  0.57  sec,  the  minimum  fault 
dimension  is  —  1.4  km  and  the  maximum  fault  dimension  is  LBUZ  —  4.1  km. 

Using  aftershock  hypocentral  locations,  Seeber  it  at.,  (1086)  suggest  a  small  source 
area  centered  at  a  depth  of  4.5  km  and  about  1  km  in  the  largest  dimension. 

Kanamori  and  Anderson,  (1075),  provided  a  theoretical  basis  for  associating  the  mt 
of  a  small  event  to  the  fault  length  (I). 


3  log  L  *  mh  -  4.0 


(14) 


From  (14)  we  obtain  a  fault  length  of  «.74  km.  The  differences  in  the  three  fault-length 
estimates  can  probably  be  attributed  to  the  assumptions  inherent  in  the  different 
methods. 

Having  defined  the  source  time  function  of  the  Ardsley  m*in  shock  with  respect  to 
its  largest  aftershock  we  now  turn  to  the  problem  of  generating  synthetic  seismograms 
for  the  event. 

A  suite  of  synthetic  seismograms  has  been  generated  for  the  Ardsley  event  at  the 
distance  and  azimuth  of  station  SRNY  for  comparison  to  the  recorded  data.  Records 
were  synthesized  at  source  depths  of  0,  2,  6,  8  and  10  km  up  to  a  frequency  of  5  Hz.  The 
New  England  velocity  model  (Taylor  it  at. ,  1880)  and  Q  structure  (Mitchell,  1881)  were 
used. 

Because  the  main  shock  appears  to  be  the  superposition  of  three  sub-events  (asper¬ 
ity  ruptures)  separated  in  time,  it  is  not  appropriate  to  compare  the  main  shock  directly 
to  the  suite  of  synthetics,  which  were  generated  ■fuming  a  simple  Green’s  function.  We 
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therefore  compere  the  aftershock  seismogram  to  the  synthetics.  To  match  the  5  Hz  high 
frequency  limit  of  the  synthetics  and  eliminate  microeeismic  noise,  the  data  were  band 
pass  filtered  between  .5  and  5  Hz  prior  to  the  comparison. 

The  vertical  component  traces  are  shown  in  Figure  30  with  the  appropriate  data 
trace  placed  at  a  depth  of  4.5  km.  The  first  three  phases  to  arrive  are  the 
Pg,  PiP ,  and  $PiP  phases  illustrated  in  Figure  40.  Phases  that  are  initially  downgoing 
such  as  PiP  arrive  earlier  as  the  source  depth  increases  and  phases  that  are  initially 
upgoing,  especially  depth  phases  such  as  tPiP,  arrive  later  as  the  source  depth  increases. 
The  match  of  the  synthetics  to  the  first  three  arrivals  of  the  data  is  quite  good  although 
the  amplitudes  are  not  well  matched.  In  addition,  approximately  3.5  secs  after  Pg  there 
is  a  clear  phase  on  the  data  trace.  Most  likely  it  is  tPmP ,  though  this  is  not  a  prom¬ 
inent  arrival  in  the  synthetics  and  synthetic  sPmP  arrives  slightly  later  for  a  source 
depth  of  4.5  km. 

Station  RSNY  at  44.548  *  N  74.530  *  W  is  400  km  from  the  epicenter  of  the  Ardsley 
event  and  at  nearly  the  same  azimuth  as  station  SRNY  (-15.4  *  vs.  -8.0  * ).  It  recorded 
both  the  main  shock  and  the  aftershock  discussed  above  but  at  a  much  lower  signal-to- 
noise  ratio.  Because  the  synthetic  match  at  station  SRNY  was  excellent  and  we  have  a 
good  model  for  the  far  field  displacement  time  history  along  this  azimuth,  we  have 
attempted  to  model  the  data  at  station  RSNY. 

A  summary  of  the  depth  estimates  for  northeastern  events  is  presented  below  in 
Table  5. 


TABLE  5 

Source  Depths  of  Northeast  U.S.  Earthquakes 


Date 

Event 

Magnitude 

Reported  Depth 

This  Study 

Oct  07, 1083 

Goodnow 

5.2 

7,14 

10 

Jan  10, 1082 

New  Hampshire 

4.5  mk 

3-11 

8 

Oct  11, 1083 

Ontario 

4.1 

12 

14 

Oct  10, 1085 

Ardsley 

4.0  m, 

5 

4.5 

Oct  02, 1085 

Quarry 

2.3  m( 

0 

0 

Aug  31, 1082 

Adirondack 

2.8  m. 

? 

7 

Oct  23,  1084 

Adirondack 

3.4  mt 

? 

12 

Oct  30, 1085 

Amsterdam 

2.7  m( 

? 

1 

May  20, 1083 

Maine 

4.2 

2 

12 

The  suite  of  synthetics  for  RSNY  were  generated  using  a  point  source  at  depths  of 
0,  2,  4,  6,  8,  10  and  12  km.  The  results  were  then  convolved  with  an  approximation  to 
the  far  field  displacement  time  history  so  that  a  comparison  could  be  made  with  the  data 
(Figure  41).  The  vertical  data  recorded  at  RSNY  in  Figure  41  are  placed  at  the  depth 
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Figure  SO.  Vertical  component  synthetic  seismograms  for  the  Ardaley,  NY,  aftershock 
and  data  at  station  SRNY,  A— 100  km. 
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determined  from  comparison  of  the  synthetics  and  data  at  SRNY  (4.5  km).  Although 
noisy,  the  fit  of  the  data  to  the  synthetics  is  good.  The  depth  phases  are  very  well 
matched  between  seconds  two  and  five  as  is  the  arrival  at  second  nine.  What  has  not 
been  modeled  is  the  high  amplitude  coda  between  these  two  sets  of  arrivals.  This  may 
be  the  effect  of  scattering,  multipathing,  or  the  simplicity  of  the  synthetic  model. 


Kuril/Kamchatka  Data 

The  Kuril/Kamchatka  data  will  be  discussed  collectively  rather  than  event  by  event 
as  was  done  for  the  northeast  U.S.  data.  The  event  and  OSS-IV  locations  are  shown  in 
Figure  7.  Sixteen  events  were  studied,  of  which  fourteen  were  assigned  depths  (Table  3). 
All  of  the  events  are  above  magnitude  4.5  and  range  in  distance  from  574  km  to  1307  km 
from  OSS-IV. 

The  OSS  instrument  was  emplaced  380  m  below  the  sea  floor  in  5467  m  deep  ocean. 
The  primary  motivations  for  placing  the  instrument  in  a  borehole  are  to  enhance  the 
coupling  of  the  instrument  and  to  remove  the  sensors  from  the  noisy  sear  bottom  environ¬ 
ment.  Soft  sediments  at  the  ocean  bottom  have  strong  shear-wave  velocity  gradients, 
which  trap  shear-wave  energy  close  to  the  sea-floor  interface.  By  installing  the  sensors 
well  below  the  sea  floor,  the  shear  wave  noise  is  avoided.  Unfortunately,  being  well 
below  the  sea  floor  interface  can  complicate  the  receiver  site  effects.  At  typical  compres¬ 
sions!  wave  sediment  velocities  (2.0  km/sec)  the  travel  time  from  the  instrument  to  the 
sea-floor  interface  and  back  is  0.38  seconds.  The  two-way  travel  time  for  the  sea-surface 
reflection  is  7.7  seconds.  The  sea-surface  reflection  is  the  strongest  of  the  two  (Figure 
42-50).  As  seen  from  the  data,  the  prominent  sea-surface  reflection  would  make  any 
depth  phase  arriving  near  7.7  seconds  difficult  to  identify.  A  pP  phase  arriving  at  7.7 
seconds  corresponds  to  an  event  depth  of  about  23  km. 

For  the  northeastern  U.S.  data,  adaptive  polarization  analysis  was  used  to  help 
identify  phases.  The  apparent  angle  of  incidence  trace  identified  arrivals  with  similar 
phase  velocities  and  the  RZ  trace  helped  identify  true  arrivals.  The  unusual  location  of 
OSS-IV,  beneath  the  sea  floor,  makes  use  of  adaptive  polarization  analysis  suspect.  A 
strong  sear  bottom  reflection  adds  downgoing  arrivals  to  the  usual  upgoing  waves.  The 
apparent  angle  of  incidence  trace  would  then  become  useless  and  the  RZ  trace  would 
indicate  both  P  and  SV  type  motion.  Event  Kl  (Figure  43  and  44)  is  an  example  of  this 
type  of  event.  During  the  initial  arrivals,  the  RZ  trace  is  both  positive  and  negative 
(usually  a  sign  of  Rayleigh  type  motion)  and  the  apparent  angle  of  incidence  trace  is 
confused.  In  contrast,  the  searsurface  reflection,  which  does  not  include  the  sea-bottom 
reflections,  has  stable  negative  apparent  angle  of  incidence  and  RZ  trace. 

Events  with  a  large  signal- to-noise  ratio,  Kl,  K3,  K4,  K5,  K6,  K7,  K10,  Kll,  and 
K18  all  exhibit  a  common  character  in  the  first  part  of  the  wavetrain.  The  first  arrival  is 
generally  small  and  is  followed  by  three  or  more  arrivals  with  increasingly  large  ampli¬ 
tudes  in  the  next  two  seconds.  The  later  arrivals  could  be  either  depth  phases  (pP,  »P ), 
source  effects,  or  receiver  site  effects.  The  reported  depths  of  the  events  eliminates  depth 
phases  as  the  source  of  these  arrivals  for  all  but  the  K7  event  (6  km  depth).  If  the 
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arrivals  are  due  to  source  effects  then  a  remarkably  consistent  source  mechanism  would 
be  required  for  all  of  these  events.  Although  source  effects  are  certainly  present  in  the 
data,  the  regularity  of  the  arrivals  over  so  many  events  strongly  suggests  that  they  are 
caused  by  site  effects.  The  additional  arrivals  are  also  present  in  the  sea-surface 
reflection. 

Uncertainty  in  the  amount  of  skew  correction,  problems  associated  with  the  location 
of  the  OSS  and  lack  of  good  long  range  synthetics  made  depth  determination  for  this 
data  set  extremely  difficult.  Four  different  seismologists  with  varying  degrees  of  observa¬ 
tional  experience  independently  evaluated  the  data  for  depth  and  obtained  consistent 
picks  for  only  one  of  the  16  events  (Table  6).  Ironically,  that  event,  K5,  has  no  pub¬ 
lished  depth.  The  depth  determined  by  the  observers  for  K5  was  about  16  km. 


TABLE  6 

Differences  in  Travel  Times  Between  the  First 

Arrival  and  the  First  Identifiable  Depth  Phase 
as  Picked  by  Four  Analysts  for  the  Kuril/Kamchatka  Data 

Event 

Analyst  #1 

Analyst  #2 

Analyst#3 

Analyst  #4 

Kl 

4.7 

2.0 

11.8 

11.0 

K3 

3.5 

3.6 

5.0 

7.4 

K4 

3.2 

5.4 

14.0 

10.0 

K5 

5/5.5 

£.4 

5.7 

5.4 

K6 

3.8 

2.0 

13.0 

3.0 

K7 

3.5 

2.0 

13.0 

13.0 

K8 

2.2 

5.4 

13.0 

1.25 

K10 

3.2/5.3 

4.6 

2/3.0 

10.7 

Kll 

4.0 

2.1/14.6 

6.0 

13.0 

K12 

5.5 

5.7 

2.1/3 

2.2 

K14 

8.0 

5.7 

.8/0/16.2 

4.1 

K15 

5.5 

4.6 

4.8 

3.0 

K16 

4.8 

4.6 

0.5/10.6 

11.3 

K18 

8.0 

5.4 

2.3/11 

6.6 

K10 

4.7 

3.2 

5.6 

0.0 

K20 

2.0 

2.0 

? 

5.2 
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Figure  42.  120  seconds  of  unfiltered  data  from  OSS-IV  for  event  Kl,  A=861  km.  Large 

bursts  on  the  East-West  trace  at  seconds  40  and  100  are  minute  marks  that 
are  used  for  skew  correction. 
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Figure  \6.  State  filtered  (sf)  and  band-pass  filtered  (bp)  P- wave  data  for  event  K4, 
A =605  km.  The  band-pass  filter  passed  frequencies  between  2  and  15  Hz. 
The  polarization  state  filter  was  applied  to  the  band-passed  data  and  passed 
rectilinear  motion. 
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Figure  48.  State  filtered  (sf)  and  band-pass  filtered  (bp)  P- wave  data  for  event  K6, 
A *948  km.  The  band-pass  filter  passed  frequencies  between  2  and  15  Hz. 
The  polarization  state  filter  was  applied  to  the  band-passed  data  and  passed 
rectilinear  motion. 


Figure  50.  State  filtered  (sf)  and  band-pass  filtered  (bp)  P-wave  data  for  event  K8, 
A —046  km.  The  band-pass  filter  passed  frequencies  between  2  and  15  Hz. 
The  polarization  state  filter  was  applied  to  the  band-passed  data  and  passed 
rectilinear  motion. 
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Figure  5J.  State  filtered  (af)  and  band-pass  filtered  (bp)  P- wave  data  for  event  K14, 
A=»753  km.  The  band-pass  filter  passed  frequencies  between  2  and  15  Hz. 
The  polarization  state  filter  was  applied  to  the  band-passed  data  and  passed 
rectilinear  motion. 
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State  filtered  (sf)  and  band-pass  filtered  (bp)  P- wave  data  for  event  K20, 
A— 1291  km.  The  band-pass  filter  passed  frequencies  between  2  and  15  Hz. 
The  polarisation  state  filter  was  applied  to  the  band-passed  data  and  passed 
rectilinear  motion. 
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SECTION  V. 

SUMMARY  AND  CONCLUSIONS 


Application  of  the  depth  phase  method  for  event  depth  determination  in  the 
northeastern  U.S.  and  the  Kuril/Kamchatka  area  gave  mixed  results.  In  the  northeast, 
excellent  results  were  obtained  (see  Table  5)  while  in  the  Kuril /Kamchatka  area  very 
poor  performance  was  recorded.  The  conditions  that  led  to  the  dramatic  contrast  in  the 
results  can  be  used  to  quantify  the  method's  strengths  and  weaknesses. 

One  major  difference  between  the  two  data  sets  is  signal-to-noiae  ratios.  With  the 
exception  of  the  Amsterdam  earthquake  and  quarry  blast,  all  of  the  recorded  signals  for 
the  northeast  events  have  good  signal- to- noise  ratios,  whereas  signal-to-noiae  ratios  for 
the  Kamchatka  data  were  poor.  This  was  due,  in  part,  to  the  difference  in  epicentral  dis¬ 
tances;  for  northeast  events  all  epicentral  distances  are  less  than  330  km,  while  for  the 
Kamchatka  data  they  are  all  between  S70  and  1301  km.  Also,  the  propagation  paths 
from  many  of  the  Kamchatka  events  travene  a  major  tectonic  boundary,  the  Kuril 
Trench,  and  this  may  contribute  to  a  decrease  in  signal- to-noise.  The  signal- to-noise 
ratio  plays  a  significant  role  in  the  depth  determinations.  Tests  of  the  polarisation  state 
filter  method  used  to  help  isolate  depth  phases  have  shown  that  performance  degrades 
significantly  below  signal- to-noise  ratios  of  1.  This  hampered  the  efforts  to  identify 
depth  phases  for  the  Kamchatka  events. 

The  second  m^jor  difference  between  the  data  sets  is  the  data  quality.  Data  col¬ 
lected  in  the  northeast  U.S.  were  recorded  by  broad  band  stations  RSNY  and  SRNY. 
Both  record  data  digitally  and  have  well  known  characteristics.  The  OSS  data,  collected 
using  a  remote  sea-bottom  recording  system,  has  problems  with  skew,  timing,  and  possi¬ 
bly  instrument  noise.  The  polarisation  methods  described  in  this  report  require  that  the 
three  components  of  motion  be  accurately  aligned  in  time.  Although  skew  was  corrected 
for  in  the  OSS  data,  some  inaccuracies  may  have  persisted. 

The  physical  location  of  the  recording  stations  also  contributed  to  differences 
between  the  data  sets.  Both  RSNY  and  SRNY  are  well  coupled  stations  at  or  near  the 
free  surface  of  the  Earth.  In  contrast,  OSS- IV  is  well  below  the  sea  floor  in  deep  ocean. 
Superimposed  upon  the  upgoing  seismic  wavefield  at  OSS- IV  are  reflections  from  the 
sear  bottom  interface  and  the  sea  surface.  The  resulting  wavefield  is  complicated,  mating 
use  of  the  adaptive  polarisation  method  for  phase  identification  of  minimal  value  as  well 
as  reducing  the  effectiveness  of  the  polarisation  state  filter. 

Another  difference  that  bears  on  the  results  is  the  way  the  two  data  sets  were 
analysed.  The  northeastern  seismograms  were  compared  to  synthetic  depth  sections, 
whereas  we  simply  attempted  to  pick  upgoing  and  downgoing  phases  on  the  complex 
Kamchatka  seismograms.  The  synthetic  seismograms  guided  our  identification  of  phases 
considerably  and,  without  them,  it  would  have  been  very  difficult  to  identify  properly- 
paired  primary  and  depth  phases  in  the  northeastern  data. 


Readout  Associates,  lac. 


Merck  1987 


AJTAC  final  Report 


The  preceding  discussion  pointed  out  possible  reasons  that  the  depth  phase  method 
gave  depths  consistent  with  reported  depths  for  the  northeastern  U.S.  data  but  did  not 
work  on  the  Kamchatka/Kuril  data.  In  the  course  of  applying  the  method  to  the 
northeastern  data,  we  found  that  several  other  conditions  were  important  for  its  success. 
The  Ardsley  event  revealed  that  source  complexity  is  possible  in  small  events  and  that 
the  phases  associated  with  the  complexity  could  be  misidentified  as  depth  phases.  Since 
identification  of  phases  can  be  aided  by  synthetic  seismograms,  it  is  important  to  have  a 
measure  of  the  validity  of  the  models.  A  models’  robustness  could  be  tested  by  success¬ 
fully  modeling  the  data  from  several  different  receiving  stations.  Unfortunately,  we  do 
not  often  have  on-scale  digital  data,  especially  three-component  data,  from  more  than 
one  station  in  the  northeast.  Nevertheless,  we  can  afford  a  certain  degree  of  confidence 
in  the  modelling  and  interpretation  of  earthquakes  we  selected  because:  1.  the  epicenters 
are  well  located  by  dense  regional  networks,  2.  some  of  the  depths  are  also  well  deter¬ 
mined  by  virtue  of  doee-in  permanent  or  temporary  stations,  3.  for  more  than  half  the 
earthquakes,  fault-plane-solutions  have  been  obtained  from  P-wave  first  motions  and  4. 
several  regional  velocity  structures  have  been  estimated  from  refraction  experiments, 
teleseismic  residuals,  and  travel-time  curves.  More  comparisons  need  to  be  made,  how-i 
ever,  to  establish  whether  the  velocity  models  so  derived  (#4  above)  are  the  most 
appropriate  for  waveform  modelling  at  the  frequencies  of  interest.  It  is  apparent  from 
the  comparisons  we  have  presented  here  that  the  data  are  not  perfectly  modelled  by  the 
synthetics.  Future  investigations  should  focus  on  improving  these  results  through  better 
understanding  of  both  data  and  synthetics.  Now  that  digital  data  are  available  from 
regional  networks  ,  there  is  great  potential  for  discovering  much  more  about  source 
mechanics  and  crustal  phase  propagation.  Future  insights  can  give  us  a  better  idea  of 
how  to  model  the  entire  source-path-receiver  system. 

In  conclusion,  the  depth  phase  method  for  depth  determination  can  work  quite  well 
under  certain  conditions.  The  problem  that  remains  to  be  solved  is  to  identify  when 
these  conditions  are  met  for  a  particular  data  set. 
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ATTNs  Dr.  L.  Turnbull 
OSI/NSD,  Room  5G48 
Washington,  D.C.  20505 

U.  S.  Arms  Control  and  Disarmament  Agency 
ATTN;  Dr.  M.  Elmer 

Verification  and  Intelligence  Bureau,  Rm  4953 
Washington,  D.C.  20451 

U.  S.  Arms  Control  and  Disarmament  Agency 
ATTNs  Mrs.  M.  Hoinkes 
Multilateral  Affairs  Bureau,  Rm  5499 
Washington,  D.C.  20451 

U.  S.  Geological  Survey 

ATTNs  Dr.  T.  Hanks 

National  Earthquake  Research  Center 

345  Middlefield  Road 

Menlo  Park,  CA  94025 


OTHER  GO  VERMONT  AGENCIES  (continued) 

0.  S.  Geological  Survey  MS-913 
ATTHt  Dr.  R.  Maaae 
Global  Seismology  Branch 
Box  25046,  Stop  967 
Denver  Federal  Center 
Denver,  CO  80225 

UNIVERSITIES 

Boston  College 
ATTN i  Dr.  A.  Kafka 
Western  Observatory 
381  Concord  Road 
Weston,  MA  02193 

California  Institute  of  Technology 
ATTHt  Dr.  D.  Harkrider 
Seisaological  Laboratory 
Pasadena,  CA  91125 

Columbia  University 
ATTHt  Dr.  L.  Sykes 

Lament -Doherty  Geological  Observatory 
Palisades,  NY  10964 

• 

Cornell  University 

ATTHt  Dr.  M.  Barazangi 

INSTOC 

Snee  Hall 

Ithaca,  NY  14853 

Harvard  University 
ATTHt  Dr.  J.  Woodhouse 
Hoffman  Laboratory 
20  Oxford  Street 
Cambridge,  MA  02138 

Massachusetts  Institute  of  Technology 
ATTHt  Dr.  S.  Solemn,  Dr.  H.  Toksos,  and  Dr.  T.  Jordan 
Department  of  Barth  and  Planetary  Sciences 
Cambridge,  MA  02139 

Southern  Methodist  University 
ATTHt  Dr.  B.  Herrin 
Geophysical  Laboratory 
Dallas,  TX  75275 

State  University  of  New  York  at  Binghamton 
ATTHt  Dr.  F.  Wu 

Department  of  Geological  Sciences 
Vestal,  NY  13901 


UNIVERSITIES  (continued) 


St.  Louis  University 

ATTN i  Or.  0.  Nuttli  and  Dr.  R.  Herrmann 
Department  of  Barth  and  Atmospheric  Sciences 
3507  Laclede 
St.  Louis  ,  MO  63156 

The  Pennsylvania  State  University 
ATTN*  Dr.  S.  Alexander 
Geosciences  Department 
403  Deike  Building 
University  Park,  PA  16802 

University  of  Arizona 
ATTN*  Dr.  T.  Wallace 
Department  of  Geosciences 
Tucson,  AZ  85721 

University  of  California,  Berkeley 
ATTN*  Dr.  T.  Me Evilly 
Department  of  Geology  and  Geophysics 
Berkeley,  CA  94720 

University  of  California  Los  Angeles 
ATTN*  Dr.  L.  Knopoff 
405  Hilgard  Avenue 
Los  Angeles,  CA  90024 

University  of  California,  San  Diego 
ATTN*  Dr.  J.  Orcutt 
Scripps  Institute  of  Oceanography 
La  Jolla,  CA  92093 

University  of  Colorado 
ATTN*  Dr.  C.  Archambeau 
CXRES 

Boulder,  CO  80309 

University  of  Illinois 
ATTN*  Dr.  S.  Grand 
Department  of  Geology 
1301  West  Green  Street 
Urbana,  IL  61801 

University  of  Michigan 

ATTN*  Dr.  T.  Lay 

Department  of  Geological  Sciences 

Ann  Arbor,  MI  48109-1063 

University  of  Nevada 
ATTN*  Dr.  K.  Priestley 
Mackay  School  of  Mines 
Reno,  NV  89557 


UNIVERSITIES  (continued) 

University  of  Southern  California 

ATTN i  Dr.  K.  Aki 

Center  for  Earth  Sciences 

University  Park 

Los  Angeles,  CA  90089-0741 

DEPARTMENT  OP  DEFENSE  CONTRACTORS 

Applied  Theory,  Inc. 

ATTN:  Dr.  J.  Trulio 
930  South  La  Brea  Avenue 
Suite  2 

Los  Angeles,  CA  90036 

Center  for  Seismic  Studies 
ATTN:  Dr.  C.  Romney  and  Mr.  R.  Perez 
1300  N.  17th  Street,  Suite  1450 
Arlington,  VA  22209 

BNSCO,  Inc. 

ATTN:  Mr.  G.  Young 
5400  Port  Royal  Road 
Springfield,  VA  22151 

ENSCO,  Inc. 

ATTN:  Dr.  R.  Kemerait 
445  Pineda  Court 
Melbourne,  FL  32940 

Gould  Inc. 

ATTN:  Mr.  R.  J.  Woodard 
Chesapeake  Instrument  Division 
6711  Baymeado  Drive 
Glen  Bumie,  MD  21061 

Pacific  Sierra  Research  Corp. 

ATTN:  Mr.  F.  Thomas 
12340  Santa  Monica  Boulevard 
Los  Angeles,  CA  90025 

Rockwell  International 
ATTN:  B.  Tittmann 
1049  Camino  Dos  Rios 
Thousand  Oaks,  CA  91360 

Rondout  Associates ,  Inc , 

ATTN:  Dr.  P.  Pomeroy 
P.O.  Box  224 
Stone  Ridge,  NY  12484 

Science  Applications,  Inc. 

ATTN:  Dr.  T.  Bache,  Jr. 

P.O.  Box  2351 
La  Jolla,  CA  92038 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (continued) 

Science  Horizons  2 

ATTN:  Dr.  T.  Cherry  and  Dr.  J.  Minster 
710  Encinitas  Blvd. 

Suite  101 

Encinitas,  CA  92024 

Sierra  Geophysics,  Inc.  2 

ATTN:  Dr.  R.  Hart  and  Dr.  G.  Mellman 
11255  Kirkland  Way 
Kirkland,  HA  98124 

SRI  International  1 

ATTN:  Dr.  A.  Florence 
333  Ravensworth  Avenue 
Menlo  Park,  CA  94025 

S -Cubed  a  Division  of  y 

Maxwell  Laboratories  Inc. 

ATTN:  Dr.  S.  Day 

P.O.  Box  1620 

La  Jolla,  CA  92038-1620 

S -Cubed,  A  Division  of  y 

Maxwell  Laboratories  Inc. 

ATTN:  Mr.  J.  Murphy 
11800  Sunrise  Valley  Drive 
Suite  1112 
Rea ton,  VA  22091 

Teledyne  Geotech  2 

ATTN:  Dr.  X.  Der  and  Mr.  W.  Rivers 
314  Montgomery  Street 
Alexandria,  VA  22314-1581 

Woodward -Clyde  Consultants  2 

ATTN:  Dr.  L.  Burdick  and  Dr.  J.  Barker 
556  El  Dorado  Street 
Pasadena,  CA  91109 


NON-US  RECIPIENTS 

Blacknest  Seismological  Center  1 

ATTN:  Mr.  Peter  Marshall 

Atomic  Weapons  Research  Establishment 

UK  Ministry  of  Defense 

Brimpton,  Reading  RG7-4RS 

United  Kingdom 

National  Defense  Research  Institute  1 

ATTN:  Dr.  Ola  Dahlman 
Stockholm  80 ,  Sweden 
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NON-OS  RECIPIENTS 


NTNP  NORSAR 

ATTN i  Dr.  Prod*  Ringdal 
P.O.  Book  SI 
N-2007  K jailor 
Norway 

OTHER  DISTRIBUTION 

To  bo  dotorainod  by  tho  projoct  offico 


